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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to 
SI  (metric)  units  as  follows: 


Multiply 
degrees  (angle) 
inches 

pounds  (force) 
pounds  (mass) 
tons  (force) 


_ 52 _ 

0.01745329 

2.54 

4.448222 

0.4535924 

8.896444 


To  Obtain 
radians 
centimetres 
newtons 
kilograms 
kilonewtons 
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BATTLEFIELD  DUST  FROM  EXPLODING  MUNITIONS:  CONTRIBUTION 


BY  CRATERING  FROM  ARTILLERY  AND  MORTAR  PROJECTILES 


PART  I:  INTRODUCTION 

Background 

1.  The  problem  of  dust  loading  in  the  near-surface  atmosphere  as  it 
affects  the  modern  battlefield  environment  was  clearly  established  during  the 
North  African  campaigns  of  World  War  II.  Both  men  and  machines  were  severely 
affected  by  it — men,  by  respiration  and  eye  irritations  and  infections,  and 
machines,  by  engine  and  abrasion  problems  that  caused  high  failure  rates  and 
compounded  the  difficult  logistics  in  that  theater. 

2.  Since  that  time,  the  level  of  technology  on  the  battlefield  has 
increased  markedly,  and  dust-laden  air  space  has  become  an  even  greater  factor 
in  military  operations,  as  demonstrated  in  the  recent  Middle  East  wars.  High- 
technology  tracking  and  aiming  devices  may  not  function  properly  during 
periods  when  critical  dust-cloud  concentrations  are  present.  Whether  the  dust 
stems  from  natural  causes,  such  as  a  dust  storm  (or  "haboob”) ,  or  from  a  vari¬ 
ety  of  battlefield  activities,  such  as  moving  vehicles  and  exploding  muni¬ 
tions,  the  problem  is  the  same  when  dust  concentrations  reach  critical  levels. 

3.  Dust  loading  resulting  from  military  operations  common  to  the  bat¬ 
tlefield  varies  widely,  depending  upon  geographic  location  (i.e.,  soil  com¬ 
position),  climate,  season,  and  to  some  extent,  even  the  time  of  day. 

Research  by  the  US  Army  Engineer  Waterways  Experiment  Station  (WES) , 

Vicksburg,  MS,  is  aimed  at  developing  realistic  methods  for  predicting  dust 
loadings  resulting  from  military  activities  in  desert  terrain,  especially  in 
those  areas  where  conventional  military  operations  appear  most  likely.  As 
input  to  a  prediction  methodology  for  determining  dust  concentrations  within 
dust  clouds  generated  by  exploding  munitions.  Science  and  Technology  Corpora¬ 
tion,  under  contract  to  WES,  compiled  a  data  base  documenting  crater  size  and 
shape  and  ejecta  volume  from  explosions  representative  of  the  range  of  energy 
yields  common  to  such  munitions.  The  surveyed  explosions  involved  bare 
(uncased)  high-explosive  (HE)  charges  as  well  as  static  and  live-fired 
munitions  in  a  variety  of  soils.  This  field  study  was  followed  by  an  analysis 
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of  the  data  base,  which  established  volume-based  soil  scaling  fac¬ 
tors  (VBSSFs)*  suitable  for  mathematical  modeling  algorithms,  as  well  as  other 
cratering  relations.  This  report  combines  the  findings  of  the  preceding 
research  efforts. 


Purposes  and  Scopes 

4.  The  purposes  of  this  study  are:  (a)  to  publish  a  comprehensive  data 
base  that  can  be  used  to  quantify  VBSSFs  as  they  vary  with  soil  type  and,  to 
the  degree  possible,  to  assess  the  effect  of  vegetative  cover  on  the  VBSSFs 
and  on  the  crater-ejection  process  for  near-surface  charges  in  the  weight 
range  of  roughly  1  to  120  kg**  (which  includes  the  range  of  explosive  weights 
common  to  most  conventional  munitions  expected  on  the  modern  battlefield)  and 
(b)  to  p  -rform  a  detailed  empirical  analysis  of  the  acquired  data  base  to 
quantify  the  VBSSFs  and  crater-ejecta  parameters  for  the  soils  and  shot  geo¬ 
metries  included. 

5.  The  third  and  final  phase  of  this  research  (Phase  III)  will,  if 
authorized,  review  and  hopefully  improve  the  numerical  input  to  various  dust 
codes,  e.g.,  the  Combined  Obscuration  Model  for  Battlefield-Induced  Contami¬ 
nants  (COMBIC)  using  the  findings  of  this  study. 


*  Crater  volume  divided  by  charge  weight  raised  to  some  exponential  value. 

**  Strictly  speaking,  the  kilogram  (1  kg  -  2.205  lb,  mass)  is  a  unit  of  mass, 
not  weight.  It  is  used  here  as  a  unit  of  force  (weight)  in  the  Earth's 
gravitational  field,  i.e.,  1  kg  =  2.205  lb,  force,  in  place  of  the  correct 

but  unfamiliar  newton  (IN*  0.225  lb,  force).  Actually  the  120-kg  weight 
limit  stated  above  represents  265  lb  or  I, 17V  N. 
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PART  II:  COMPILATION  PROCEDURES 

Literature  Survey 

6.  Resources  consulted  and  documented  in  the  accomplishment  of  this 
study  are  listed  in  the  References  at  the  end  of  the  main  text.  This  listing, 
which  is  alphabetical  by  principal  author  (or  organization) ,  includes  material 
that  furnished  direct  input  to  this  study.  Other  material  used  for  general 
background  information  but  not  cited  in  the  text  are  listed  in  the 
Bibliography. 


Data  Tabulation  Formatting 

7.  Narrative  descriptions  of  most  of  the  test  sites  from  which  data  for 
this  study  were  used  are  provided  in  Part  III.  The  results  of  the  literature 
search  are  discussed  in  Part  IV,  and  the  nine  data  tables  found  at  the  end  of 
the  report  contain  VBSSFs  for  all  explosion-crater  events  considered  as  well 
as  the  results  of  crater-ejecta  measurements,  calculations,  and  estimates. 

The  events  recorded  in  these  tables  are  mostly  uncased  charges,  but  Tables  5 
and  9  include  data  from  static  and  live  projectile  firings.  Finally,  there 
are  five  "nondata"  tables,  four  of  which  document  the  development  of  study 
parameters  used  in  connection  with  the  projectile  craters  and  one  of  which 
records  the  development  of  crater-ejecta  relations  for  uncased  charges. 

8.  Upon  completion,  the  data  base  matrix  of  Tables  1-9  (row  by  column) 
was  516  by  35;  each  row  is  identified  by  a  unique  number  for  ease  of  rtier- 
ence.  The  data  tables  are  contained  in  both  the  Phase  I*  and  Phase  II** 
reports,  as  well  as  here,  with  final  additions  and  corrections  being  included 
in  the  Phase  II  report  and  here.  However,  there  are  a  number  of  incomplete 
entries,  limited  by  the  measurements  that  could  be  found  or  reasonably 
inferred. 


*  Allen  D.  Rooke,  Jr.,  and  John  N.  Strange.  1986.  "Battlefield  Dust  from 
Exploding  Munitions:  A  Data  Base  Descriptive  of  the  Contribution  by 
Cratering,"  STC  Technical  Report  2096,  Science  and  Technology  Corporation, 
Hampton,  VA. 

**  John  N.  Strange,  and  Allen  D.  Rooke,  Jr.  1986.  "Battlefield  Dust  from 
Exploding  Munitions:  An  Analysis  of  the  Phase  I  Data  Base,"  STC  Technical 
Report  2123,  Science  and  Technology  Corporation,  Hampton,  VA. 
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PART  III:  SOIL  DESCRIPTIONS  OF  THE  VARIOUS  TEST  AREAS 


9.  One  of  the  major  objectives  of  this  study  was  to  develop  VBSSFs 
that  could  be  correlated  directly  to  specific  soil  types,  such  as  those 
specified  by  the  Unified  Soil  Classification  System  (USCS).  Unfortunately,  it 
was  found  during  the  course  of  the  literature  search  that  only  in  the  more 
recent  studies  of  cratering  had  soil  descriptions  for  various  test  sites  been 
sufficiently  examined/tested  to  state  a  USCS  designation.  Where  such  designa¬ 
tions  were  made,  obviously  they  were  used.  Where  designations  were  not  made 
within  the  body  of  a  given  report,  the  authors  of  this  report  have  estimated 
USCS  designations  based  on  word  descriptions  of  the  natural  soils  present  at 
the  various  sites  and,  in  some  cases,  on  personal  knowledge  of  soil  conditions 
at  various  sites.  Estimated  USCS  designations  have  been  clearly  identified  as 
such. 

10.  For  the  ourpose  of  this  study,  where  so  many  site  descriptions  are 
nonspecific,  a  triangular  classification  chart  (Figure  1),  similar  to  that  of 
the  US  Department  of  Agriculture  (USDA) ,  has  been  prepared  to  give  a 
quantitative  indication  of  the  proportion  of  sand,  silt,  and  clay  present  in 
the  natural  soils.  Note  that  the  soil  types  identified  in  Figure  1  are 
uniquely  related  to  the  word  descriptions  given  for  the  various  sites. 

!1.  In  addition,  a  delineation  of  soil  type  relative  to  particle  size 
is  given  below: 


Material  Particlp  Size  (PS),  mm 


Clay 

PS 

< 

0.005 

Silt 

0.005 

<  PS 

< 

0.07 

Sand  (fine) 

0.07 

<  PS 

< 

0.20 

(medium) 

0.20 

<  PS 

< 

0.6 

(coarse) 

0.6 

<  PS 

£ 

2 

Gravel  (fine) 

2 

<  PS 

< 

10 

(coarse) 

10 

<  PS 

< 

70 

Boulders 

PS 

> 

70 

These  particle-size  statistics  are  based  on  a  concensus  evaluation  of  a  number 
of  different  size  groupings  from  different  organizations,  specif ically,  the 
American  Society  for  Testing  and  Materials,  the  Federal  Aviation  Administra¬ 
tion,  the  USCS,  The  American  Association  of  State  Highway  Officials,  the 
USDA,  and  the  Massachusetts  Institute  of  Technology. 
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12.  The  word  descriptions  of  the  soils  at  the  various  sites,  with  their 
stated  or  estimated  USCS  designations,  are  presented  in  the  following  para¬ 
graphs,  along  with  references.  Where  soil  descriptions  or  even  test  site 
locations  were  nonspecific,  the  authors  have  substituted  their  best  judgments. 
Descriptions  of  vegetative  cover  are  also  included  where  applicable  and/or 
available.  In  most  cases  where  there  is  no  description  of  vegetative  cover, 
the  reader  may  assume  that  testing  was  done  on  bare  ground  or  that  vegetation 
was  so  sparse  that  it  could  be  ignored.  However,  descriptions  for  the  follow¬ 
ing  sites  were  not  available:  Fort  Benning  (paragraph  21),  Panama  Canal 
(paragraph  37),  and  miscellaneous  sites  (paragraph  41). 
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Big  Black  Test  Site  (BBTS) 

Cratering  effects  of  surface  and 
buried  HE  charges  in  loess  and  clay 

13.  The  BBTS  is  a  WES  test  site,  located  about  12  miles  east-southeast 
of  Vicksburg,  Miss.  The  loess  test  pad  at  this  site  had  an  average  moisture 
content  of  19  percent  during  the  test  period.  Other  properties  were:  average 
density  of  1,640  kg/m"^  (-102  lb/ft^),  liquid  limit  of  44.8,  a  plastic  limit  of 
24.3,  and  a  plasticity  index  of  20.5.  The  USCS  designation  was  estimated  as 
CL,  and  the  percent  sand/silt/clay/moisture  content/degree  of  satura¬ 
tion  (SSCMS)  is  estimated  as  10e/40e/50e/10e/20e  (the  lower-case  "e"  is  used 
here  to  indicate  estimated  quantities) . 

14.  The  clay  test  pad  had  an  average  moisture  content  of  21  percent 
during  the  test  period.  Other  mechanical  properties  were:  average  density  of 
1,874  kg/m^  (-117  lb/ft'*),  liquid  limit  of  43,  a  plastic  limit  of  23,  and  a 
plasticity  index  of  20.  The  estimated  USCS  designation  was  CL-CH;  the  SSCMS 
is  5e/20e/75e/10e/25e  (WES  1958a). 

Craters  formed  by 

small  explosions  in  dry  sand 

15.  The  test  bed  for  the  series  was  a  dry  (moisture  content  less  than 
f.  percent),  well-graded  sand.  The  USCS  designation  was  SW;  the  SSCMS  is 
l00e/0e/Ge/<6/10e  (Sager  1962). 

Effects  of  a  soil- 

rock  interface  on  cratering 

16.  These  tests  were  conducted  in  a  uniformly  graded  sand  overlying  a 
massive  concrete  block  of  sufficient  thickness  to  act  as  solid  rock  for  the 
scale  of  the  experiment.  Only  those  cratering  data  from  sand-layer  thick¬ 
nesses  of  sufficient  depth  to  negate  the  eifects  of  the  underlying  concrete 
block  are  used  here. 

17.  Moisture  contents  within  the  sand  layer  ranged  from  4.4  percent 
near  the  surface  to  7.5  percent  at  the  sand-concrete  interface;  it  was  clas¬ 
sified  as  dry-to-moist ,  with  an  estimated  USCS  designation  of  SP.  The  SSCMS 
is  100e/0e/0e/4e  to  8/15e.  For  the  wet-sand  tests,  the  moisture  content  was 
12  to  15  percent  with  saturation  of  about  20  percent  (WES  1958b). 
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The  influence  of  a 

shallow  water  table  on  cratering 

18.  A  special  test  bed  of  uniform  washed  masonry  sand  was  used  for 
these  tests.  Only  the  data  from  shots  where  the  water  table  was  sufficiently 
deep  so  as  not  to  affect  crater  formation  are  used  here.  The  USCS  designation 
was  SP,  with  an  SSCMS  of  100/0/0/15e/25e  (Carnes  1981). 

Eglin  Air  Force  Base 

19.  Two  separate  series  of  explosion  tests,  SMOKE  WEEK  II  and  SMOKE 
WEEK  III,  were  conducted  at  Test  Range  C5A  in  the  eastern  part  of  Eglin  Air 
Force  Base.  The  soil  at  this  site  was  classified  as  a  more  or  less  uniform 
sand  with  minor  amounts  of  silt  and  clay.  The  USCS  designation  was  estimated 
to  be  SP,  while  the  SSCMS  designation  is  80e/10e/10e/<9/20e .  The  moisture 
content  at  depths  of  3  to  10  cm  averaged  in  the  range  of  5  to  7  percent. 

20.  Vegetation  consisting  of  grasses,  herbs,  and  forbs  covered  about 
50  percent  of  the  area.  Vegetation  was  removed  from  shot  sites  for  SMOKE 
WEEK  II,  but  was  not  removed  for  SMOKE  WEEK  III  (Mason  and  Long  1983a). 

Fort  Benning,  GA 

21.  The  Portable  Bunker  Tests  were  conducted  on  the  Coolidge  Mortar 
Range  and  the  Spafford  Artillery  Range  at  Fort  Benning.  The  soil  at  both 
ranges  was  a  moist  silty  sand  with  a  USCS  designation  of  SM.  On  the  mortar 
range,  a  clayey  sand  (USCS  designation  SC)  was  found  at  a  depth  of  0.6;  how¬ 
ever,  this  horizon  was  penetrated  only  once  by  the  explosion  trials.  The 
SSCMS  designation  is  60e/30e/10e/5e/7e  (Hoot  1971). 

Fort  Carson,  CO 


22.  The  surface  soil  at  Site  1  of  the  Dust  Obscuration  Test  Series  I 

(DOT  I)  was  a  sandy  clay  material,  light  brown  in  color.  The  USCS  designation 

was  a  sandy  clay  (CL).  Near-surface  moisture  contents  taken  in  the  upper  10 
to  15  cm  ranged  from  6.7  to  26.7  percent,  the  higher  value  occurring  soon 
after  a  light  rain.  Cone  index  readings  (a  measure  of  the  bearing  capacity  of 

soils  relating  mainly  to  trafficability)  at  depth  were:  50  at  surface,  200  at 
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15  cm,  and  750+  at  30  cm.  The  SSCMS  designation  is  35e/10e/55e/7  to  25/20e  to 
50e. 

23.  At  Site  2  (DOT  II)  the  soil  was  a  sandy  clay,  sandy  silt  classified 
as  CL.  At  the  time  of  the  testing,  the  moisture  content  ranged  from  7.2  to 
17.5  percent,  the  higher  value  again  occurring  soon  after  a  light  rain.  Cone 
index  readings  at  depth  were:  100  at  the  surface  and  750+  at  15  to  30  cm. 

For  Site  2,  the  SSCMS  designation  is  30e/35e/35e/7  to  17.5/15e  to  30e. 

24.  Some  ground  cover  was  present:  "short  grasses,  blue  grama,  Russian 
thistle,  lamb's  quarters,  and  prickly  pear."  Ground  coverage  was  50  to 

70  percent  on  the  average;  plant  height  generally  ranged  from  10  to  20  cm 
(Long,  Mason,  and  Durst  1984;  Long  and  Williams  1985;  Long  et  al.  1985). 

Fort  Polk,  LA 


Munitions /bare  charge 
equivalence  Series  II  (MBCE  II) 

25.  These  experiments  were  conducted  within  the  confines  of  Range  37  at 
Fort  Polk,  as  were  the  Battlefield  Environments  in  Tailored  Soils  (BETS) 
tests.  The  soils  to  a  depth  of  2  m  consisted  of  0.3  m  of  tan  sandy  clay  over- 
lying  about  1  m  of  tan  to  reddish-tan  plastic  clay,  which  in  turn  overlay 
reddish-tan  to  brown  clay.  The  dominant  soil  type,  within  the  depths  affected 
by  the  various  shots,  was  classed  as  a  sandy  clay  (estimated  USCS,  CL),  with 
an  SSCMS  of  about  30e/10e/60e/10  to  30  (average  ~20)/20e  to  60e  (average  40e) 
(Mason  and  Long  1983b) . 

BETS 

26.  The  BETS  tests  were  conducted  at  a  time  when  the  average  moisture 
content  was  lower  than  on  the  MBCE  program,  or  about  10  percent.  Especially 
prepared  test  beds  were  used  for  these  tests.  They  were:  (a)  inorganic 

clay  (CH) ;  (b)  mixture  of  inorganic  silts,  some  clay,  and  some  fine  sand  (ML); 
(c)  poorly  graded  sand  (SP);  (d)  inorganic  clay-sand  mixture  (CH/SP); 

(e)  silty  sand  (ML-SP) ;  and  (f)  kaolin  (CH) . 

27.  The  organic  content  of  the  Range  37  noils  varied  from  1+  to  nearly 
7  percent,  with  an  average  of  about  3.5  percent  (Mason  and  Long  1983b). 
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Mono  Lake ,  CA 


28.  A  number  of  cratering  experiments  were  conducted  on  the  south  shore 
of  Mono  Lake,  where  the  soil  was  pumice  sand,  with  only  minor  amounts  of  silt 
and  clay.  The  sand  was  uniform  in  size,  with  grain  size  in  the  main  a  mil¬ 
limetre  or  less  in  diameter.  The  USCS  designation  was  SP,  with  an  estimated 
SSCMS  of  90e/5e/5e/varies/varies  (Davis  1967). 

Tropic  Test  Center  (TTC) ,  Isthmus  of  Panama 

29.  Explosion  tests  were  conducted  at  three  sites  in  TTC  under  the  Bat¬ 
tlefield  Obscuration  in  the  Tropics  test  program.  At  TTC  Empire  Range  6  on 
the  Pacific  side  of  the  Isthmus,  the  dominant  soil  was  a  moist  silty  clay  with 
a  liquid  limit  that  ranged  from  50  to  60.  The  USCS  designation  was  nearly 
always  MH;  the  SSCMS  averages  about  24e/55e/21e/35e/70e. 

30.  The  Mindi  Farm  test  area  is  on  the  Atlantic  side  of  the  Isthmus; 
the  soil  there  was  primarily  silt  with  some  fine  sand  and  a  trace  of  clay. 

This  material  was  moist-to-wet  (but  generally  more  wet  than  moist);  moisture 
contents  ranged  from  70  percent  at  the  surface  to  40  percent  at  a  depth  equal 
to  crater  depth.  The  USCS  designation  was  MH,  with  an  estimated  SSCMS  that 
averages  about  38e/57e/5e/60/85e. 

31.  At  the  Mindi  Farm  site,  charges  were  detonated  in  three  different 
levels  of  vegetation:  bare  soil  cleared  of  all  vegetation,  Gynevium 
sagittatum  0.3  to  0.5  m  high,  and  the  same  vegetation  3  to  4  m  high. 

32.  The  Pina  Beach  site  is  located  on  the  Atlantic  side  of  the  Isthmus, 
approximately  1  km  southwest  of  the  mouth  of  the  Chagres  River.  The  beach 
material  was  fine  coastal  sands  with  traces  of  silt  and  gravel.  Charges  were 
detonated  in  three  soils  and  three  types  of  vegetation.  The  soils  were  "white 
saturated  sand  (shoreline);  white,  wet  (top  centimetre  partially  dry)  sand; 
black,  wet  (top  centimetre  partially  dry)  sand;  the  types  of  vegetation  were 
Ipomoea  pes-aaprae  (morning  glory),  Hymenoaallis  americana  (spider  lily),  and 
Paniaum  maximum  (2  to  3  m  high."  The  average  SSCMS  for  Pina  Beach  is  esti¬ 
mated  as  93e/7e/0/10e/20e  (Martinucci  and  Fuchs  1981). 
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White  Sands  Missile  Range  (WSMR) 


Dusty  Infrared  Test 
(DIRT)  Series  I  and  II 

33.  The  DIRT  I  tests  were  conducted  in  the  Orogrande  (southeast)  area 
of  WSMR.  The  near-surface  soils  in  the  area  were  composed  principally  of 
drv-to-moist  brown  silty  sand.  The  estimated  USCS  was  classified  as  SM;  the 
SSCMS  is  taken  as  60e/35e/5e/3  to  13/20e. 

34.  The  DIRT  II  series  was  conducted  in  the  Queen  15  area  in  north- 
central  WSMR.  The  material  was  an  alluvial  deposit  described  as  a  moist  silty 
clay,  estimated  to  be  CL-ML,  with  silt  and  some  sand.  The  water  table  lay  at 
a  depth  of  about  2  m.  The  SSCMS  is  15e/35e/50e/10  (near  the  surface;  at  1-m 
depth  the  moisture  content  was  24  percent,  for  an  overall  average  10  to 
15e)/20e  to  40e  (average  *  30e). 

Battlefield  Induced 
Contaminants  (BICT)  Series  III 

35.  This  series  was  conducted  in  the  Orogrande  area,  WSMR,  described  in 
paragraph  33. 

Mun itions/bare 

charge  equivalence  Series  I 

36.  Series  I  of  MBCE  was  held  in  the  Queen  15  area  of  WSMR,  described 
in  paragraph  34.  (MBCE  II  is  discussed  in  paragraph  25.)  (See  Joachim  and 
Davis  >1Q83);  Rooke  (1983);  Williams  and  Lebron  (in  preparation)). 

Miscellaneous  Sites 


Panama  Canal  studies 

37.  The  Panama  Canal  tests  were  conducted  in  residual  clay  estimated  to 
be  OL  (wet)  and  in  marine  muck  estimated  to  be  OH,  near  saturation.  Moisture 
contents  in  both  materials  were  in  excess  of  40  percent.  Specific  weights  for 
the  residual  clay  were:  1,767  kg/m^  (wet)  and  1,169  kg/m^  ('dry);  the  SSCMS  is 
5e/20e/75e/40  to  50e/70e  to  80e  (US  Special  Engineering  Division,  The  Panama 
Canal  1948a  and  1948b). 
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Small  explosion 
tests — Project  MOLE 

38.  These  experiments  were  conducted  In  four  rather  basic  soil  types: 
dry  clay  (Utah),  sand  with  some  gravel  (Nevada  Test  Site),  wet  sand  (Cali¬ 
fornia),  and  moist  clay  (California).  The  USCS  classifications  were  estimated 
as  follows  (Sachs  and  Swift  1955): 


Word  Description 
Dry  clay 

Sand  with  some  gravel 
Wet  sand 
Moist  clay 

Suffield  Experimental 
Station  (SES) ,  Alberta,  Canada 

39.  The  Distant  Plain  Series  was  conducted  within  the  confines  of  the 

Watching  Hill  and  Drowning  Ford  Test  Ranges,  SES,  where  the  general  topography 

was  that  of  a  gently  rolling  plain  sparsely  vegetated  by  buffalo  and  spear 

grasses.  The  soil  at  the  Watching  Hill  site  was  a  glacial  deposit  with  strata 

as  defined  by  the  following  tabulation: 


Estimated 

USCS 

Designation 

CL 

GM 

SW 

CL 


Word  Description _ 

Silty  clay 

Lean  clay,  sandy  or  silty 
Clay 

Coarse  sand  and  gravel 


Estimated 

USCS 

Designation 

OL 

SM 

CL 

GP 


Depth,  m 
0  -  1.5 

1.5  -  9 

9  -  25 

25  -  ? 


The  SSCMS  is  10e/25e/65e/10e/15e  through  the  0-  to  3-m  horizon  (Jones, 
Spackman,  and  Winfield  1959;  Jones  et  al.  1970). 

AO.  The  origin  of  the  soils  at  the  Drowning  Ford  test  area  was  also  a 
glacial  deposit,  with  the  underlying  strata  described  as  follows: 
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Word  Description 
Clayey  sil t 
°andy  s:.lt 
‘'and/  »ravei  mix 
Cl nyey  silt 


Estimated 

uses 

Designation 

ML 

SM 

SW 

ML 


Depth,  m 
0  -  3.5 

3.5  -  4.0 

4.0  -  13.5 
3.5  -  ? 


For  the  Drowning  Ford  site,  the  SSCKS  as  a  rough  average  is  10e/55e/35e/ 
10e/15e  within  the  0-  to  3-m  horizon  (Jones,  Spackman,  and  Winfield  1959; 

Jones  et  al.  1970). 

Effects  of  underground  explosions 

41.  These  tests  were  conducted  at  three  sites:  Princeton,  N.  J.;  Clear 
Lake,  Tex.;  and  Natchez,  Miss.  The  estimated  USCS  designations  for  the  three 
sites  were  as  follows:  Princeton — MH,  underlain  at  shallow  depth  by  CL;  Clear 
Lake — CL;  and  Natchez — CL.  There  are  insufficient  data  avail,  hie  on  these 
tests,  conducted  in  the  1950s,  to  permit  an  estimate  of  SSCMS  (Engineering 
Research  Associates  1952). 

Underground  explosion  test  program 

42.  Tests  were  conducted  in  rather  commonly  occurring  soil  types, 
namely,  dry  clay  (estimated  as  CL),  dry  sand  (estimated  as  SP-SM  dry),  and  wet 
clay  (estimated  as  CL).  Here  again,  no  SSCMS  can  be  estimated  (Engineering 
Research  Associates  1952). 

Operation  Prairie  Flat 

43.  This  test  was  conducted  at  SES  Watching  Hill  Range,  described  in 
paragraph  39  (Jones,  Spackman,  and  Winfield  1959;  Jones  et  al.  1970). 
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PART  IV:  ANALYSIS 


Purposes 

44.  The  main  efforts  in  the  analysis  of  the  compiled  data  are  to  docu¬ 
ment  in  a  quantitative  fashion  the  VBSSFs,  the  scaling  law  relating  crater 
volume  and  the  quantity  of  explosive  (charge  weight) ,  and  to  develop  means  of 
estimating  ejecta  volume  as  a  fraction  of  crater  volume.  In  so  doing,  differ¬ 
ences  between  cased  and  uncased  charges  have  been  noted.  The  analysis  of  the 
compiled  data  is  reflected  in  the  remaining  portions  of  this  report. 

Determination  of  Crater  Volume:  Crater  Shape  Factors 


Apparent  craters 

45.  Figure  2  is  an  Idealized  profile  of  a  crater  formed  by  a  near¬ 
surface  explosion.  If  one  assumes  that  the  apparent  crater  volume  is  directly 
proportional  to  a  variation  of  the  square  of  the  crater  radius  times  the 
crater  depth,  then  it  is  possible  to  develop  a  crater  shape  factor  in  order  to 
calculate  crater  volume  from  the  measures  of  radius  and  depth.  The  simplest 
way  to  arrive  at  such  a  factor  is  to  compare  the  actual  apparent  crater  volume 
with  that  of  a  fictitious  right  circular  cylinder  having  the  same  radius  and 

height  (in  this  case  depth)  as  the  crater,  i.e.,  compare  the  apparent  crater 

2 

volume  V  with  irr  d  ,  where  r  is  the  apparent  crater  radius  and  d  is 
a  a  a  a  a 

the  apparent  crater  depth. 

2 

46.  If  Vfl  is  assumed  to  be  directly  proportional  to  irrflda  »  an 
assumption  borne  out  by  experimental  results,  then 


V 

a 


(1) 


and  for  the  true  crater,  denoted  by  the  subscript  t  , 


vt  -  V-tV 


(2) 
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Figure  2.  Idealized  half-crater  profiles 


By  definition,  the  C  terms  become  the  crater  shape  factors  K  used  in 

s 

Tables  1-9  for  both  apparent  and  true  craters.  In  the  tables,  apparent  and 
true  shape  factors  are  separated  by  a  slash,  apparent/ true.  In  the  text  dis¬ 
cussion,  they  will  be  shown  as  (K  )  and  (K  ) 

s  a  s  t 

47.  The  data  from  fully  reported,  bare-charge  cratering  events  in  var¬ 
ious  soils  from  various  heights  (depths)  of  burst,  and  with  charge  weights 
ranging  from  0.5  kg  to  as  much  as  500,000  kg  (TNT  equivalent),  resulted  in  a 
(K  )  that  averaged  0.45  with  a  standard  deviation  of  o  of  only  0.07.  The 

S  3 

range  of  (K  )  values  was  from  0.28  to  0.59.  With  (K  )  «*  0.45  ,  the  crater 

s  a  s  a 

shape  more  nearly  approximates  a  paraboloid  (Kg  *  1/2)  as  opposed  to  a  cone 
(K  *  1/3)  or  a  hyperboloid  (K  ■  2/3).  Even  though  the  crater  aspect  ratio 


(r  / d  )  for  the  data  base  developed  varied  from  1.5  to  10,  the  crater  shape 

o  O. 

factor  varied  only  from  0.28  to  0.59,  with  the  mean  value  being  0.45. 

48.  With  the  crater  shape  factor  known,  it  is  possible  to  calculate  the 

apparent  crater  volume  from  the  measured  crater  depth  and  radius.  With  (K  ) 

s  a 

-0.45  ,  the  apparent  crater  volume  is  given  by 

V  =  (K  )  7ir2d  =  0.45  (rr2d  )  (3) 

a  s  a  a  a  V  a  a/  v  y 


or 


V  =  1.41  r2d 
a  a  a 

Thus  for  bare-charge  apparent  craters,  a  shape  factor  of  0.45  should  be  used 
in  calculating  crater  volume  from  measured  values  of  radius  and  depth  for  the 
range  of  charge  positions  considered,  which  may  be  expressed  in  terms  of  an 
equivalent  spherical  TNT  charge  with  radius  ,  where  -2.5  <  <  4.0  . 

True  craters 

49.  A  similar  treatment  of  true  crater  results  shows  (K  )  to  have  a 

s  t 

value  of  0.40,  indicating  a  tendency  to  be  a  bit  more  conical  in  shape  than 

the  apparent  crater.  A  typical  true  crater  shape  is  shown  in  Figure  2  along 

with  the  apparent  crater  and  the  elastic-  and  plastic-response  zones.  The 

compression/compaction  cavity  (the  bottom  portion  of  the  true  crater  profile) 

is  obviously  paraboloidal,  while  the  side  walls  (upper  portion  of  the  true 

crater  profile)  form  a  truncated  cone,  thus  making  the  overall  shape  of  the 

true  crater  a  combination  of  the  two,  i.e.,  a  cone  (K  =  1/3)  and  a  paraboloid 

s 

(K  =  1/2). 

s 

50.  With  (K  )  =  0.40  ,  and  knowing  the  true  crater  radius  (r  )  and 

St  t 

depth  (d^),  it  is  possible  to  calculate  the  true  crater  volume  (V^)  from 

vt  -  <Ks)t("rtdt)-°-4oKdt)  <4> 


or 


V 

t 


1.26  r^dt 
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Thus,  for  true  craters,  a  shape  factor  of  0.40  should  be  used  to  calculate 
true  crater  volume  from  measured  values  of  radius  and  depth,  again  within  the 
range  of  charge  positions  considered  (-2.5  ^  <  4.0). 

Crater  shape  factors,  munitions 

51.  Implicit  in  this  study  is  a  requirement  to  compare  the  limited  data 
base  available  for  exploding  munitions  with  that  for  bare  charges.  The  pur¬ 
pose,  of  course,  is  to  answer  the  question:  How  well  do  craters  formed  by 
bare  charges  simulate  those  formed  by  commonly  used  battlefield  munitions, 
such  as  artillery  and  mortar  projectiles?  Bare  charges  are  safer  and  more 
convenient  than  munitions  for  conducting  battlefield  cratering  (and  ejecta/ 
dust-lofting)  experiments,  but  only  insofar  as  they  satisfactorily  represent 
the  cratering  effects  of  munitions. 

52.  A  typical  live-fire  apparent  crater  profile  for  an  artillery  pro¬ 
jectile  launched  by  cannon  is  shown  in  Figure  3.  Craters  from  live-fire 
munitions  differ  little  from  static-fire  craters,  where  shells  are  manually 
emplaced,  even  though  in  the  live-  ire  case,  shell  fragments  possess  a  forward 
velocity  vector  equal  to  the  shell's  velocity  vector.  Shell  craters  are 
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typically  asymmetric,  with  the  uprange  radius  generally  less  than  the  down- 
range  radius;  both  of  these  differ  from  the  lateral  (cross-range)  radii,  which 
are  generally  the  same . 

53.  Figure  4  shows  apparent  crater  shape  factors  for  munitions.  In 
spite  of  the  asymmetry  of  munition  craters  compared  with  the  usually  symmet¬ 
rical  bare-charge  craters,  the  shape  factors  for  munition  craters  range  from 
0.22  to  0.78  as  read  from  the  abscissa  values,  excluding  the  seemingly 
anomalous  point  P  ;  this  gives  an  average  value  of  K  =  0.46  .  Thus,  in  the 


9  ai-uu  uortar 
<3  4.2-IN  UORT AR 
■  I  05  MU  HOW  |  T  2  £  R 
•  I5S-MU  HOWITZER 


O  61-MU  MORTAR 
®  4,2  *  I  N  MORTAR 

O  lOi-uw  HOWITZER 
O  I  55- MM  HOWITZER 
A  122-UM  GuN-HOWIT  Z.E  R 
V  I  52  -MM  CUN-HOWITZER 


LIVE  EIRE  O 
-J 

o 

1  i 


23 


RIGHT  CIRCULAR  CYLINDER 


absence  of  a  uniquely  defined  munition  shape  factor,  the  average  radius  and 
measured  depth  may  be  used  in  conjunction  with  the  appropriate  bare-charge 
shape  factor  (0.45  for  apparent  and  0.40  for  true)  to  calculate  crater  volume. 


VBSSFs 


Apparent  craters: 
unrestricted  analysis 

54.  Dimensional  analysis  considerations  dictate  that  crater  volume, 
both  apparent  and  true,  should  vary  directly  as  charge  weight,  l.e.,  a  dou¬ 
bling  of  charge  weight  should  result  in  the  doubling  of  the  crater  volume  for 
similarly  scaled  charge  positions  and  similar  soils.  Thompson  and  DeVore 
(1982),  using  cratering  data  from  a  number  of  sources,  found  that,  for  charge 
yields  appropriate  to  their  study,  the  apparent  crater  volume  varied  as  the 
charge  weight  to  the  1.111  power,  i.e.. 


V  « 
a 


W 


1.111 


(5) 


where  W  is  charge  weight. 

55.  To  test  the  validity  of  Equation  5  to  the  data  base  compiled  during 
the  Phase  I  portion  of  this  study,  a  plot  was  made  in  which  restricted  (W 
<  1,000  kg)  and  unrestricted  (W  <  1,000,000  kg)  graphs  were  drawn  (Figure  5). 
The  plot  is  an  all-data  plot,  meaning  that  data  from  all  different  soil  types 
and  all  charge  positions  were  included.  The  plot  exhibits  a  vertical  scatter 
of  roughly  two  orders  of  magnitude,  implying  that  the  extremes  of  soil  types 
(from  wet  clay  to  rock)  and  the  extremes  of  charge  position  (-2.3  <  R£  <  4) 
are  capable  of  producing  only  a  two-order-of-magnitude  scatter  in  apparent 
crater  volume  for  a  given  charge  weight.  Thus,  the  effect  of  wide  variations 
in  soil  types  as  well  as  variations  in  charge  position  in  the  near-surface 
regime  may  be  expected  to  produce  little  more  than  a  two-order-of-magnitude 
scatter  in  plots  of  apparent  crater  volume  versus  charge  weight  over  a  wide 
range  of  yields. 
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56.  The  graph  of  Figure  5,  in  which  approximately  95  percent  of  the 

478  data  points  fall  within  the  scatter  band  shown,  provides  a  massive  array 
of  points  for  slope  determination.  A  least  squares  fit  of  the  data  suggests 
that 


W 


0.85 


(6) 


When  charge  yields  are  restricted  such  that  W  s  1,000  kg  ,  the  variation  of 
apparent  crater  volume  with  charge  weight  follows  the  proportional  equation 


V  a  W 
a 


0.84 


(7) 


In  both  Equations  6  and  7,  the  exponent  is,  for  all  practical  purposes,  the 
same  and  is  substantially  below  the  1.111  variation  reported  by  Thompson  and 
DeVore  (1982)  from  their  less  extensive  data  base. 

57.  Comparison  of  Equations  6  and  7  convincingly  shows  that  even  to  the 
million-kilogram  yield  level  the  same  scaling  laws  are  applicable  as  when 
charge  weights  are  restricted  to  the  thousand-kilogram  level.  This  similarity 
indicates  that  there  is  no  significant  change  in  cratering  phenomena  or  in 
energy  partitioning  phenomena  over  a  seven-order-of-magnitude  change  in  charge 
yield. 

58.  Figures  6  and  7  present  plots  of  apparent  crater  volume  as  a  func¬ 
tion  of  charge  weight  for  0  <  R  <1  and  -1.0  <  R£  <  0  ,  respectively.  The 
plots  show  slopes  (exponents  of  W  )  of  0.84  and  0.90  when  W  <  1,000,000  kg 
and  0.78  and  0.89  when  W  5  1,000  kg.  Again,  the  vertical  scatter  is  roughly 
two  orders  of  magnitude  for  both  plots. 

59.  It  appears  at  this  point,  by  virtue  of  the  preponderance  of  data 
shown  in  Figures  5-7,  that  the  1.111  exponent  previously  assumed  is  too  high. 
These  data  tend  to  support  an  exponent  in  the  range  of  0.85  to  0.90  for  charge 
weights  less  than  1,000,000  kg. 

60.  It  must  be  emphasized  at  this  point  that  the  exponents  derived  from 
Figures  5-7  are  for  correlations  where  all  soil  types  and  all  charge  positions 
or  ranges  of  charge  positions  are  lumped  together  to  provide  the  respective 
plots.  This  lumping  of  data  presumes  that  each  soil  type  and  each  restricted 
range  of  charge  position,  as  an  overall  average,  will  follow  the  same  scaling 
laws  and  that  the  plots  for  more  restrictive  correlations  will  be  described 
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Variation  of  apparent  crater  volume  with  charge  weight  (all  soils  types 
“1*0  5  Rr  5  0).  Graph  contains  137  data  points 


only  by  different  ordinate  intercepts,  i.e.,  by  the  multiplying  constant  C 
in  the  basic  equation 


V  -  cwn 
a 


(8) 


Restricted  analysis 

61.  As  a  practical  matter,  the  most  likely  occurring  charge  position 
for  a  munition  is  one  near  the  surface,  or  air-ground  interface,  i.e.,  where 
R  “  0  .  Figure  8  shows  a  plot  of  apparent  crater  volume  versus  charge  weight 
for  charges  positioned  such  that  -0.5  £  Rc  £  0.5  .  The  plot  shows 


V  * 
a 


W 


0.96 


(9) 


Figure  9  shows  a  similar  but  more  restrictive  plot  regarding  charge  position; 

here,  the  center  of  gravity  of  the  charge  is  at  the  surface  (R  =0).  This 

c 

relationship  is  described  by  the  proportional  equation 

V  <*  W0,94  (10) 

a 

The  vertical  scatter  band  in  Figure  8  is  roughly  two  orders  of  magnitude, 
which  is  principally  a  measure  of  the  effect  of  soil-type  variations  as  well 
as  small  variations  in  charge  position.  In  Figure  9,  where  the  charge  posi¬ 
tion  is  fixed,  the  vertical  scatter  is  just  slightly  greater  than  one  order  of 
magnitude,  inferring  that  the  extremities  of  soil  types  cause  roughly  an 
order-of-magnitude  variation  in  apparent  crater  volume  when  the  charge  posi¬ 
tion  is  held  constant. 

62.  All  apparent  crater  volume  plots  thus  far  (Figures  5-9)  demonstrate 

that 


V  «  Wn 
a 


(ID 


where 


0.78  <  n  ^  0.96 
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MUNITION  DATA 
BARE-CHARGE  DATA 


W-Charge  Weight,  kg  (TNT  equivalent) 

Variation  of  apparent  crater  volume  with  charge  weight  (all  soil  types 
-0.5  <  R  <  0.5).  Graph  contains  113  data  points 


MUNITION  DATA 
BARE-CHARGE  DATA 


Variation  of  apparent  crater  volume  with  charge  weight  (all  soil  types 
R  *  0) .  Graph  contains  90  data  points 


The  average  value  of  n  ,  giving  equal  weight  to  each  plot,  Is  0.88.  A 
weighted  average,  based  on  the  number  of  plotted  points  in  each  figure,  is 
0.86.  Considering  the  scatter  in  the  data  at  hand  and  the  paucity  of  data  for 
specific  values  of  R£  ,  it  seems  appropriate  to  calculate  VBSSFs  for  apparent 
crater  volume  using  0.87,  or  an  average  of  0.88  and  0.86  as  the  exponent  of 
charge  weight  rather  than  the  1.111  that  is  currently  used.  Here  n  is 
stated  to  only  two  significant  figures,  considered  the  maximum  warranted  by 
the  variability  within  the  data  base. 

Effect  of  charge 
position:  similar  soils 

63.  The  effect  of  varying  the  charge  position  while  holding  soil  type 
relatively  constant  is  shown  in  Figures  10-13,  where  VBSSFs  for  apparent 
crater  volume  are  plotted  as  a  function  of  charge  position  for  sands  (Fig¬ 
ure  10),  for  cohesive  soils  (Figure  11),  for  silty  soils  (Figure  12),  and  for 

various  rocks  (Figure  13).  Note  that  in  Figures  10-13  the  VBSSFs  (ordinate 

3  0  87 

scale)  are  presented  in  units  of  m  /W  rather  than  the  previously  used 

exponent  on  W  of  1.111.  The  dashed  portion  of  each  curve,  representing 
values  of  J  2.5  ,  is  a  "best  trend"  guess. 

64.  If  the  upper  bounds  of  the  plots  presented  In  Figures  10-13  are 
extrapolated  to  a  value  of  R£  -  -10  ,  or  the  burst  depth  that  maximizes 
apparent  crater  volume,  then  the  apparent  crater  volumes  that  would  be 
obtained  for  the  various  soil  types  are  given  in  the  tabulation  below  for  a 
100-kg  charge  and  are  compared  with  the  stated  volume  peaks  given  by  Thompson 
and  DeVore  (1982)  for  a  charge  of  equal  weight. 


Type* 

Sand 

Clay 

Silt 

Rock 


Calculated  Apparent 
Crater  Volume,  m 
(from  Figures  10-13) 

55 

219 

87 

5.5 


Thompson  and  Devore  (1982) 


V*1-111 

Stated  Peak  Value** 

0.373 

1.68 

0.553 

0.042 


Calculated  Apparent 
3 

Crater  Volume,  m 

62.2 

280 

92.2 

7.0 


*  Sand  *  dry  to  slightly  moist;  clay  =  wet;  silt  =  moist  to  near  wet;  rock 
=  medium  hardness.  .. 

**  Values  are  in  in  /kg 

t  Obtained  by  averaging  values  for  wet  sand  and  moist  cohesive  soils  and 
dry-to-moist  sandy  soils. 


65.  At  best,  the  comparison  discussed  in  paragraph  63  and  exhibited  in 
the  above  tabulation  is  highly  approximate.  In  the  first  place,  the  raw  data 
from  which  Thompson  and  DeVore  (1982)  determined  their  peak  values  are  not 
available  for  direct  comparison  by  plotting.  Second,  the  soil  types,  although 
qualitatively  similar,  may  not  be  sufficiently  similar  to  warrant  direct  com¬ 
parison.  Nevertheless,  the  figures  obtained  differ  by  less  than  30  percent. 

True  craters: 
unrestricted  analysis 

66.  The  "all  data"  plot  (no  delineation  as  to  soil  type  or  charge  posi¬ 
tion)  of  true  crater  volume  as  a  function  of  charge  weight  is  shown  in  Fig¬ 
ure  14.  As  seen  from  the  least  squares  fit  equations,  the  slope  of  the  "all 
weights"  graph  is  1.02.  For  this  graph,  the  true  crater  volume  for  the  DIAL 

3 

PACK  Event  (V  ■  16,300  m  ;  W  =  454,000  kg)  was  included.  When  DIAL  PACK  is 
omitted  and  the  charge  weight  is  thus  restricted  such  that  W  s  1,000  kg  ,  the 
slope  is  somewhat  steeper,  with  a  value  of  1.09.  In  Figure  14,  the  scatter  is 
less  than  two  orders  of  magnitude  for  a  given  charge  weight. 

Restricted  analysis 

67.  Figures  15-18  present  plots  similar  to  Figure  14,  but  for  specific 
ranges  of  charge  position: 


Figure 


Range  of  Charge  Position 


0  <  R  <1 
c 

R  =  0 
c 

-0.5  <  R  <0.5 
c 

-0.1  <  R  <0 
c 


In  these  figures,  the  slopes  average  1.03.  On  a  weighted  basis  (again  based 
on  the  number  of  plotted  points  in  each  graph),  the  average  is  also  1.03.  In 
these  figures,  the  scatter  is  approximately  one  order  of  magnitude  for  a  given 
charge  weight. 

68.  Over  the  range  of  charge  weights  of  chief  concern  to  this  study  (W 
S  1,000  kg),  the  proportional  equation  is  applicable: 


V  cc  w 
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Variation  of  true  crater  volume  with  charge  weight  (all  soil  types 
all  charge  positions) .  Graph  contains  69  data  points 


BARE-CHARGE  DATA 
ALL  WEIGHTS 
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Variation  of  true  crater  volume  with  charge  weight  (all  soil  types 
0  <  s  1).  Graph  contains  42  data  points 


MUNITION  DATA 
BARE-CHARGE  DATA 


Variation  of  true  crater  volume  with  charge  weight  (all  soil  types 
R  =0).  Graph  contains  21  data  points 


where 


0.96  jn  <1.09 

The  average  value  of  n  ,  as  determined  from  Figures  14-18  on  a  weighted 
basis,  is  1.03.  As  with  the  apparent  crater  data,  it  is  felt  that  the  quoting 
of  a  value  on  n  to  more  than  two  significant  numbers  is  not  warranted.  Thus 
VBSSFs  for  true  craters  should  be  determined  on  the  basis  of  n  =  1.0  ,  i.e., 

Vt 

(VBSSF)  =  (13) 

c  Wn 

where  n  is  assigned  the  value  of  unity.  This  means  that  for  true  craters, 
volume  is  directly  proportional  to  charge  weight,  i.e.,  a  doubling  of  charge 
weight  will  have  the  effect  of  doubling  the  crater  volume. 

69.  Figure  19  presents  a  plot  of  the  true  crater  VBSSF  as  a  function  of 
Rc  without  regard  to  soil  type.  Lack  of  true  crater  data  precluded  the 
development  of  plots  for  specific  soil  types  as  was  done  for  apparent  craters. 
Although  not  well  defined,  there  is  a  tendency  for  moist-to-wet  clay  data  to 
be  in  the  upper  portion  of  the  plot,  followed  by  silts  and  sands,  with  rock 
values  plotting  near  the  bottom  portion  of  the  data  envelope. 

Effect  of  Vegetation  on  Cratering 

70.  Quantification  of  the  effect  of  vegetative  cover  on  the  size  of  the 
apparent  crater  is  possible  only  in  the  case  where  cratering  experiments  were 
conducted  at  the  TTC  and  for  a  limited  number  of  experiments  conducted  by  WES. 
Description  of  vegetative  coverage  at  other  test  sites  was  not  adequate  to 
permit  comparisons  with  any  degree  of  confidence. 

71.  Craters  formed  by  explosions  on  bare  ground  as  compared  with  like 
explosions  on  or  in  the  3ame  soil  type  but  with  grass  cover  of  height  less 
than  30  cm  are  generally  about  10  percent  larger  by  volume.  Craters  formed  on 
bare  ground  (as  compared  with  those  formed  in  ground  with  vegetation  con¬ 
sisting  of  taller  grasses,  bushes,  undergrowth,  and  small  trees,  with 
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coverage  approaching  100  percent)  are  in  the  range  of  15  to  30  percent  larger 
by  volume  for  similar  soil  types. 

72.  Ejecta  ranges  are  decreased  where  concentrated  taller  vegetation  is 
present;  however,  the  ejecta  areal  density  at  close-in  ranges  is  increased. 

Effect  of  Soil  Moisture  on  Crater  Volume 

I 

73.  Unfortunately,  there  are  relatively  few  cratering  experiments  in 
which  soil  moisture  content  is  adequately  described.  This  fact,  combined  with 
data  scatter,  prevents  correlation  of  changes  in  crater  size  with  specific 

|  moisture  contents.  Qualitatively  though,  it  was  found  that  soils  that  were 

wet  but  not  saturated  to  the  point  that  slope  failures  occurred  around  the 
crater  wall  generally  produced  craters  slightly  larger  than  dry  or  moist 
soils.  For  sand,  the  linear  dimensions  of  a  crater  are  about  5  percent  larger 
for  wet  than  for  dry  or  moist  conditions.  For  silts  and  clays,  the  linear 
dimensions  of  craters  are  typically  about  10  percent  larger  in  wet  soils  than 
in  dry  to  moist  soils.  These  increases  in  linear  dimensions  translate  into 
roughly  a  15-percent  increase  in  crater  volume  for  craters  in  sand  when  wet 
and  dry-to-moist  volumes  are  compared  and  roughly  25-  to  30-percent  increase 
in  volume  when  a  similar  comparison  is  made  for  clays  and  silts. 

74.  Cratering  in  rock,  at  the  scale  of  interest  in  this  study,  is  not 
affected  by  moisture.  Shock  transmission  in  rocks  in  which  pore  spaces  are 
filled  with  water  is  significantly  increased,  but  available  data  do  not  permit 
an  evaluation  of  how  the  improvement  in  shock  transmission  affects  crater 
size.  At  the  yield  levels  of  interest,  it  is  doubtful  that  craters  in 
so-called  moist  rock  would  be  any  different  from  those  in  dry  rock. 

Comparison  of  True  and  Apparent  Crater  Volumes 

75.  The  manner  in  which  the  ratio  of  true  crater  volume  to  apparent 

crater  volume  varies  with  charge  position  is  reflected  in  Figure  20.  Although 

the  plot  shows  considerable  scatter,  the  true  crater  volume  approaches  the 

apparent  crater  volume  as  R£  becomes  increasingly  positive  and  most  likely 

converges  at  an  R  value  near  5  or  6. 
c 
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Munitions  Craters 


76.  Crater  shape  factors  for  munitions  were  discussed  in  paragraphs  52 
and  53,  where  it  was  shown  that,  despite  the  asymmetry  typical  of  munition 
craters,  there  is  little  overall  difference  between  munitions  and  bare-charge 
craters  when  crater  volumes  are  used  as  the  basis  of  comparison.  Figures  5-9 
show  that  apparent  crater  volumes  from  exploding  munitions  fall  within  the 
scatter  band  of  bare  charges.  However,  closer  inspection  reveals  some  dif¬ 
ferences  between  the  two,  as  well  as  explicit  differences  associated  with  the 
various  munitions  and  their  detonation  geometries.  Some  of  the  differences 
noted  may  be  due  to  differences  in  the  soils,  especially  moisture  content. 
Unfortunately,  as  with  most  crater  data,  the  scatter  tends  to  mask  these  dif¬ 
ferences  and  any  statistical  significance  that  might  be  associated  with  them. 
The  following  paragraphs  give  a  detailed  discussion  of  the  munitions  portion 
of  the  data  base  compiled  in  this  study. 

Comparison  of  Crater  Volumes  for  Munitions  and  Bare  Charges 

77.  Figure  21  compares  volumes  for  nonmunition  (bare-charge)  data  with 
those  for  munitions  only.  The  bottom  graph  (Figure  21b)  is  taken  from 
Tables  5  and  9,  with  a  total  of  146  data  points  representing  195  shots  (note 
that  one  point  in  Table  9  averages  50  shots).  Six  munitions  are  represented, 
each  showing  considerable  vertical  scatter.  Part  of  this  scatter  is  due  to 
the  range  of  detonation  geometries.  While  the  '’live’’  fire  of  Table  9  is  all 
by  "fuze  quick,"  i.e.,  no  built-in  delay,  the  static  detonations  of  Table  5 
include  several  different  projectile  orientations,  both  above  and  below 
ground.  Figure  21b  may  be  thought  of  as  representing  a  band  of  cratering 
results  from  artillery  and  mortars  fired  for  various  effects  ranging  from 
antipersonnel  (fuze  quick)  to  destruction  of  fortifications  (fuze  delay). 

78.  Perhaps  the  most  severe  restriction  concerning  the  munitions  data 
base  is  evident  in  Figure  21b,  and  that  is  the  narrow  domain  of  TNT-equivalent 
charge  weights  along  the  x-axis.  This  restriction,  combined  with  the  y-axis 
scatter,  reduces  confidence  in  the  computed  data  fit.  In  view  of  the  small 
variation  in  x-values,  four  types  of  curve  fits  were  examined  in  detail; 

Table  11  lists  equations  and  associated  statistics  for  these  curves.  Of  the 
four,  the  fit  Y  ■  A  +  BX  consistently  gave  the  highest  Index  of 
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Determination  (I)  (see  Table  11  for  explanation;  briefly,  1=1  is  the  best 
fit  possible)  in  this  and  subsequent  munitions  graphs.  Visual  examination  of 
this  fit,  however,  shows  that  it  is  not  satisfactory  in  the  lower  range  of 
charge  weights.  Similarly,  two  other  candidate  curves  were  rejected  as  poor 
fits  at  the  extremes  of  the  data  and  thus  not  suited  for  extrapolation  to 
higher  or  lower  yields.  With  these  considerations  and  limitations,  the  curve 

g 

Y  »  AX  was  chosen  as  the  best  fit  overall  for  this  and  subsequent  graphs, 
despite  its  lower  I-value,  It  should  be  noted  here  that  all  values  of  I  for 
these  fits  are  low,  due  to  data  scatter;  as  can  be  seen  in  Table  11,  I 
ranged  from  0.105  to  0.473. 

79.  In  Figure  21b,  the  results  obtained  with  two  munitions  bear  further 

comment:  these  are  the  Soviet  122-mm  projectile  and  the  4.2-in.*  mortar 
projectile.  The  122-mm  shell,  fired  statically  in  several  geometries,  gave 
larger  craters  than  would  have  been  predicted  from  the  data  fit.  The  reason 
for  this  is  not  known,  but  may  have  to  do  with  the  fragmentation  character¬ 
istics  of  the  shell.  The  4.2-in.  projectile  provided  craters  grouped  into  two 
distinct  sizes,  those  well  above  the  line  from  a  surface-tangent-below  geom¬ 
etry  (buried  one  projectile  length)  and  those  far  below  the  line  from  live  or 
simulated  live  (fuze  quick)  firings,  where  2  ^  <  4  .  The  craters  for  the 

two  conditions  differ  roughly  by  a  factor  of  30,  which  is  mainly  attributable 
to  differences  in  shot  geometry.  Although  the  simulated  live  firings  would  be 
expected  to  have  smaller  craters,  they  appear  unusually  small  when  compared 

g 

with  the  81-mm  projectile.  These  departures  from  the  Y  =  AX  fit  should  be 
borne  in  mind  when  predicting  craters  for  these  munitions. 

80.  In  most  cases,  live-fire  data  should  give  the  best  representation 
of  cratering  that  would  be  expected  on  the  battlefield.  However,  the  data 
base  here  is  small  (34  points  representing  83  shots),  and,  of  course,  shot-to- 
shot  control  is  inferior  to  static  firings.  The  live  firings  of  Table  9  are 
graphed  in  Figure  22a,  with  statistics  included  in  Table  12.  Note  that  the 
U.S.S.R.  122-  and  152-mm  projectiles,  which  could  only  be  fired  statically, 
have  not  been  included.  In  an  attempt  to  improve  upon  these  limitations, 
static  firings  of  all  munitions  that  appear  to  adequately  simulate  live  fire 
were  added  from  Table  5,  bringing  the  total  to  106  data  points  representing 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  7. 
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155  shots;  this  combination  of  live  plus  selected  static  firings  is  shown  in 
Figure  22b.  Excluded  are  the  events  in  Table  5  with  entries  in  the  MO  D/A 
column  of  S,  STB,  1/3BUR90,  and  1/3BUR45.* 

81.  Artillery  projectiles  are  of  particular  interest  here,  both  because 
of  their  widespread  use  on  the  battlefield  and  because  sufficient  data  are 
available  to  permit  a  closer  look.  In  Table  12,  105-  and  155-mm  howitzer 
projectiles  fired  live  or  in  a  geometry  that  simulated  live  fire  were  sub¬ 
jected  to  additional  scrutiny.  Several  comments  are  appropriate. 

a.  In  comparing  live  and  static  firings  at  the  same  test  sites, 
the  155-mm  static  firings  appeared  to  produce  larger  craters 
than  the  155-mm  live  firings,  but  results  were  mixed  for  the 
105-mm  live  and  static  firings. 

b.  In  the  case  of  the  105-mm,  a  projectile  orientation  represent¬ 
ing  the  terminal  ballistic  path  for  a  longer  range,  say,  20  deg 
as  opposed  to  10  deg,  may  result  in  a  slightly  larger  crater. 
For  the  155-mm,  the  opposite  seems  to  be  true.  It  is  specu¬ 
lated  that  these  differences  are  due  to  fragmentation  patterns. 

c.  The  more  saturated  soils  appear  to  produce  larger  craters,  as 
would  be  expected;  however,  this  cannot  be  quantitatively 
established  with  the  data  available. 

d.  The  155-mm  shell  is  roughly  2.75  times  as  heavy  as  the  105-mm, 
both  in  terms  of  total  weight  and  charge  weight.  However,  its 
apparent  crater  volume  for  live  or  simulated  live  fire  is  only 
about  1.7  times  as  large  as  that  of  the  105-mm. 

g 

82.  In  Figure  23,  the  Y  =  AX  curve  fits  from  the  preceding  dis¬ 
cussion  (paragraphs  77-80)  are  compared.  As  can  be  seen,  the  munitions  curves 
are  consistently  steeper  in  slope  than  the  "all-data"  fit,  but  are  not  sig¬ 
nificantly  displaced  from  it  in  the  charge-weight  range  of  this  study.  In 
view  of  the  limited  charge-weight  range  for  munitions,  as  well  as  theoretical 
considerations  that  seem  to  preclude  a  slope  greater  than  unity,  it  appears 
that  the  all-data  curve  provides  a  satisfactory  estimate  of  munitions-crater 
volume.  A  word  of  caution  should  be  added  regarding  extrapolation  of  muni¬ 
tions  data,  especially  upward:  the  data  base  of  this  study  consists  of  rela¬ 
tively  heavily  cased  projectiles;  the  prediction  methodology  may  not  serve  as 
well  to  predict  craters  from  thinly  cased  munitions,  such  as  general-purpose 
bombs . 


*  See  key  to  column  headings  for  Tables  1-9  for  explanation. 
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bare-charge  data  fits.  The  "all  data"  and  "all  data  less  munitions" 
graphs  are,  for  practical  purposes,  coincident 


Ejecta  Quantities 

83.  The  second  major  objective  of  this  study  was  to  survey  all  avail¬ 
able  data  on  the  cratering  effects  of  point-detonating  artillery  and  mortar 
shells  and  on  uncased  (bare)  charges  that  reasonably  simulate  the  effects  of 
these  shells,  with  the  purpose  of  tabulating  those  parameters  which  relate 
explosive  charge  size,  burst  geometry,  and  cratered  medium  to  the  quantity  of 
earth  material  ejected  from  the  crater.  An  assessment  of  this  ejecta  should 
hopefully  lead  to  a  better  estimate  of  the  fraction  of  this  material  that  com¬ 
prises  the  dust  cloud  resulting  from  such  explosions.  Since  there  are 
distinct  differences  between  the  detonations  of  HE  projectiles  and  bare 
charges,  tests  involving  the  projectiles  themselves  should  provide  the  best 
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data  available,  and  this  section  examines  the  ejecta  data  available  for  pro¬ 
jectiles.  Furthermore,  live-fired  projectiles  and  those  detonated  statically 
in  a  geometry  intended  to  simulate  live  fire  are  considered  separately,  since 
some  differences  in  cratering  between  these  two  cases  have  been  tentatively 
identified  (Figure  23).  As  will  be  seen,  only  one  live-fire  event  is  suitable 
for  analysis  of  ejecta. 

84.  Despite  volumetric  similarities,  there  are  basic  differences  in  the 
cratering  mechanisms  of  cased  and  uncased  charges  and,  to  a  lesser  extent, 
between  live  and  static  firings  of  projectiles.  Since  these  differences 
affect  the  ejection  mechanism,  they  are  briefly  discussed  here. 

a.  Bare-charge  craters  are  formed  principally  by  the  explosive 
shock  wave,  which  tends  to  crush  and  compact  the  cratered 
medium,  and  by  the  expanding  explosion-gas  "bubble,"  which  fol¬ 
lows  immediately  behind  the  shock  wave,  scouring  loosened  mate¬ 
rial  and,  in  the  case  of  a  buried  or  partially  buried  charge, 
heaving  earth  material  from  the  crater.  By  contrast,  a  cased- 
charge  crater  is  mostly  the  result  of  the  kinetic  energy  of  the 
case  fragments  as  they  dislodge  the  soil.  It  has  been  shown 
that  over  80  percent  of  the  charge  energy  is  used  in  rupturing 
the  casing,  and  this  neglects  the  strain  energy  that  is 
expended  as  the  casing  expands  to  its  limit  before  rupture 
(Batchelor  1963) .  The  fact  that  craters  from  bare  and  cased 
charges  are  comparable  in  size  for  similar  charge  weights  and 
positions  indicates  that  the  kinetic  energy  of  the  fragments  is 
roughly  equal  to  the  energy  "envelope"  of  the  bare  charge  in 
the  cratering  process. 

b.  The  asymmetries  typical  of  projectile  craters  have  already  been 
discussed  (see  Figure  3).  Reasons  for  these  asymmetries  are 
fairly  obvious  when  consideration  is  given  to  the  position  of 
the  explosive  charge  at  the  time  of  detonation,  the  projectile 
orientation,  and  the  fragmentation  pattern.  Apart  from  the 
differences  in  profiles,  craters  from  bare  and  cased  charges 
differ  in  three  respects:  (1)  lack  of  a  well-defined,  symmet¬ 
rical  lip  for  the  cased  charge;  (2)  the  fact  that  the  bare- 
charge  crater  displays  significant  plastic  flowage  as  evidenced 
by  its  upthrusted  lip  (Figure  2) ,  while  the  cased  charge  does 
not;  and  (3)  the  fact  that  fallback  material  is  not  as  notice¬ 
able  in  the  cased-charge  crater.  The  loosened  material  is 
there,  however,  as  a  thin  covering  of  the  floor  and  walls  of 
the  crater.  The  overall  effect  is  to  complicate  the  task  of 
assessing  ejecta  quantities  for  projectile  craters. 

c.  As  to  craters  from  live  and  statically  fired  projectiles,  there 
is  one  major  difference:  the  live-fired  shell  with  fuze-quick 
setting  imparts  to  its  fragments  a  velocity  component  along  the 
projectile  path,  which  does  not  occur  in  the  static  firing.  As 
the  live-fired  shell  impacts  the  ground,  fuze  detonation  occurs 
within  a  matter  of  5  to  10  ysec,  during  which  time  the  fuze  is 
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penetrating  the  earth.  Although  deceleration  has  begun,  the 
projectile  retains  most  of  its  terminal  velocity  at  the  time  of 
detonation.  For  an  artillery  shell,  this  may  be  on  the  order 
of  300  to  400  m/sec;  for  a  mortar  shell,  about  half  this  veloc¬ 
ity.  The  forward  velocity  component  thus  added  to  the  shell 
fragments  probably  enhances  true-crater  size,  when  compared 
with  a  static  projectile  placed  so  as  to  duplicate  the  live 
shell  at  the  instant  of  detonation.  This  increase  may  not  be 
noticeable  in  the  apparent  crater,  however;  overall  differences 
in  crater  sizes  between  the  two  conditions  are  usually  small. 

85.  It  is  worth  noting  here  that,  for  both  live-  and  static-fire  cases, 
many  fragments  impact  beyond  the  confines  of  the  crater,  and  these  fragments 
are  an  additional  source  of  dust  not  included  in  a  crater  study.  The  effect 
of  such  fragments  may  be  observed  in  an  artillery  low  airburst  over  dry  soil, 
where  no  cratering  occurs  but  where  a  large  dust  cloud  is  formed. 

Live-fire  projectile  data 

86.  Table  9  contains  data  on  live-fire  tests  known  to  the  authors. 
Included  are  data  from  81-mm  and  4.2-in.  mortars  and  105-  and  155-m  howitzers. 
While  the  table  is  intended  to  be  mostly  self-explanatory,  some  comments  on 
its  composition  are  in  order.  Most  importantly,  it  should  be  noted  that 
nowhere  in  any  of  the  data  for  live  or  static  firings  were  crater-ejecta 
weights  or  volumes  explicitly  given.  However,  there  are  a  number  of  tests  in 
which  sufficient  data  were  furnished  to  permit  calculation  or  at  least  rea¬ 
sonable  inference  of  ejecta  weight  or  volume.  These  calculations  or  infer¬ 
ences  then  became  the  basis  for  certain  parametric  relations  that  were 
extended  to  other  test  results.  In  the  live-fire  table,  the  experimental 
investigation  of  Shot  T-27  provided  enough  information  for  a  volumetric  bal¬ 
ance  of  the  crater  (explained  in  the  following  paragraph),  the  only  such  bal¬ 
ance  that  could  be  accomplished  for  live  fire. 

87.  If  crater-ejecta  volumes  are  not  measured  directly,  as  by  sampling, 
they  may  be  inferred  by  "balancing."  Crater  volumes  may  be  balanced  by  the 
following  expression: 


Vt  " 
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where 


V  =  true  crater  volume  (see  Figure  3) 

V  *  volume  of  ejecta 

Vfb  ”  volume  of  fallback  material 
Vc  *  volume  created  by  compaction 

=  volume  created  by  plastic  flowage  into  the  crater  lip 
As  stated  here,  all  volumes  are  in  situ  (in  the  undisturbed  state).  Fallback 
is  usually  measured  in  its  bulked  form,  and  experience  has  shown  that  reduc¬ 
tion  of  this  as-measured  volume  by  a  factor  of  about  1.5  approximates  its 
in  situ  volume.  The  grouping  (V  +  )  indicates  that  these  volumes  are 

seldom  measured  separately.  Equation  14  permitted  a  balance  of  the  volumes  of 

Shot  T-27,  from  which  the  relation  V  =  0.22  V  was  derived,  V  being  the 

e  a  a 

apparent  crater  volume.  As  will  be  shown,  this  relation  is  probably  too 
small. 

88.  The  explosive  filler  used  in  a  projectile  may  differ  significantly 
from  the  usual  energy  standard,  TNT,  and  thus  should  be  considered  in  any 
assessment  of  c ratering  results.  The  projectile  explosive  charge  in  Table  9 
is  listed  by  actual  type  explosive  and  weight  and  by  TNT  equivalence.  In  this 
case,  TNT  equivalence  is  taken  as  that  for  blast  and  impulse,  since  it  is 
these  properties  that  contribute  to  fragmentation  of  the  shell  case.  For 
example,  the  81-mm  mortar  projectile  contains  0.95  kg  of  Composition  B,  which 
is  1.10  times  as  effective  as  TNT  in  blast  and  impulse  and  is  therefore  con¬ 
sidered  equal  to  1.05-kg  TNT  for  this  study  (Headquarters,  Department  of  the 
Army  1955)  . 

89.  To  visualize  the  geometry  of  the  explosive  charge  in  relation  to 

ground  surface,  an  equivalent  charge  radius  R£  based  on  a  sphere  of  TNT  of 

equivalent  yield  is  used,  as  it  was  in  the  discussion  of  bare  charges.  This 

facilitates  comparison  of  shot  geometries.  Actually,  projectile  charges  are 

roughly  cylindrical  in  shape,  but  the  transformation  to  a  sphere  in  concept 

does  not  appreciably  affect  shot  geometry.  Using  an  average  specific  weight 

for  TNT  (about  1,522  kg/m3  or  95  lb/ft3),  R£  =  0.0539  W1/3  ,  where  W 

(charge  weight)  is  in  kilograms  and  R£  is  in  metres.  For  the  small  charges 

involved  in  this  portion  of  the  study,  a  more  practical  approach  might  be  to 

1/3 

consider  R  in  centimetres,  in  which  case  R  =  5.393  W  .  Thus,  height 
c  c 

(depth)  to  the  center  of  gravity  (CG)  of  the  cased  explosive  charge  can  be 
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expressed  in  terms  suitable  for  comparison  with  bare  charges  for  the  various 
projectile  orientations.  Actual  CG  of  the  charge  within  each  projectile  is 
estimated  from  scaled  drawings  (Headquarters,  Department  of  the  Army  1967). 

90.  Shape  factors  used  in  calculation  of  projectile-crater  volumes 
(and,  thus,  inference  of  ejecta  volumes)  were  discussed  in  paragraph  53  and 
shown  in  Figure  4.  Unlike  the  bare-charge  craters,  these  shape  factors 
exhibit  a  wide  variation  for  different  munitions,  sites,  and  shot  geometries 
and  for  live  and  static  firings. 

Statically  fired  projectiles 

91.  As  mentioned  earlier,  live-fire  tests  are  somewhat  difficult  to 
control  and  measure,  and  static  firings  are  often  substituted  to  obtain 
crater/cloud  measurements.  For  this  purpose,  it  is  necessary  to  assume  a  pro¬ 
jectile  geometry  that  would  occur  in  live  fire.  The  angle  of  impact  can  be 
taken  from  firing  tables  and  is  found  to  be  around  15  to  20  deg  for  ranges 
(distances)  at  which  howitzers  frequently  fire  and  70  to  80  deg  for  mortars. 

If  fuze-quick  action  is  assumed,  then  penetration  of  the  fuze  into  the  earth 
can  be  estimated  as  less  than  1  cm  at  the  time  of  detonation.  This  means  that 
the  body  of  the  howitzer  projectile  (at  the  hase  of  the  fuze)  is  very  nearly 
in  contact  with  the  ground  when  detonation  occurs  and  that  the  mortar  shell  is 
several  centimetres  aboveground.  For  the  howitzer  shell,  this  convenient 
assumption  means  that  the  fuze  can  be  removed,  a  blasting  cap  and  booster 
charge  inserted,  and  the  shell  placed  on  the  ground  and  tilted  to  the  desired 
angle  for  the  static  test. 

92.  In  general,  the  comments  applicable  to  live-fired  projectiles  apply 
also  to  static  firings.  There  are,  however,  some  differences  and  some  addi¬ 
tional  considerations.  Table  5  contains  static-fire  data  on  the  same  (US) 
projectiles  as  Table  9  and,  in  addition,  data  on  Soviet  122-  and  152-mm  pro¬ 
jectiles.  Information  on  the  U.S.S.R.  projectiles  was  not  as  complete  as  for 
the  US  shells,  and  calculations  of  charge  position  may  be  less  precise.  The 
range  of  test  geometries  included  is  from  simulated  fuze  quick  (referred  to  as 
"surface  tangent")  to  one  in  which  the  projectile  is  buried  to  a  depth  equal 
to  its  own  length  along  an  assumed  trajectory  path,  or  "surface  tangent  below" 
ground.  Inclusion  of  buried  or  partially  buried  projectile  tests  provides 
cratering  data  for  cases  in  which  fuze  functioning  is  slower  than  the  "super- 
quick"  action  of  modern  US  fuzes.  A  range  of  assumed  impact  angles  is  also 


included.  For  the  WSMR  tests,  where  no  projectile  shape  factor  could  be 


correlated,  the  bare-charge  shape  factor  discussed  in  paragraphs  45-48  was 
used  to  calculate  crater  volume. 

93.  Probably  the  most  Important  difference  between  Tables  5  and  9  lies 
in  the  assumption  of  the  ve/^a  relation  for  the  estimation  of  ejecta  vol¬ 
umes.  This  fraction  was  calculated  wherever  possible  for  both  munitions  and 
bare-charge  craters  and  was  used  to  arrive  at  ejecta  volumes  where  no  better 
method  could  be  found;  it  is  shown  in  Tables  1-9  under  "Ejec  Cal."  An  unusu¬ 
ally  complete  experiment  of  the  US  Army  TTC  enabled  both  direct  calculation  of 
ejecta  volumes  (explained  in  paragraph  94)  and  volumetric  balancing  on  a  num¬ 
ber  of  static-fire  munitions  craters.  Table  13  lists  important  parameters  and 

statistics  resulting  from  these  calculations.  Similar  V  /V  relations  were 

e  a 

calculated  for  bare  charges  and  will  be  discussed  later. 

94.  The  direct  calculation  of  ejecta  volumes  mentioned  previously  is 
made  possible  by  the  sampling  of  ejecta  deposited  around  some  of  the  craters 
in  the  TTC  tests.  From  these  samples,  average  "areal  density"  6  was  cal¬ 
culated  for  each  sampling  range  R  ,  excluding  organic  material  in  the  sam¬ 
ples.  By  means  of  a  least  squares  fit,  areal  density  for  the  entire  ejecta 
field  could  be  expressed  in  the  form 

6  -  C  R_m  (15) 

where  C  is  a  constant  and  the  exponent  m  is  the  slope  of  the  logarithmic 
straight-line  fit.  Equation  15  is  probably  the  most  used  approach  to  quanti¬ 
fication  of  ejecta  from  samples.  It  can  be  easily  integrated  to  solve  for 
ejecta  weight  Wg  as  follows: 


W 

e 


R-m)  R  dR 


(16) 


where  the  limits  a  ,  b  are  the  limits  of  the  ejecta  field.  The  crater  rim 

(r  )  is  usually  substituted  for  a  ,  while  the  outer  limit,  if  not  known  (as 
d 

is  the  case  here),  can  be  taken  as  infinity  with  little  error.  The  TTC  report 

included  soil  dry  densities,  from  which  wet  densities  were  calculated  and  used 

to  convert  W  to  V 
e  e 
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95.  Unfortunately,  a  problem  that  was  encountered  In  the  use  of  Equa¬ 
tions  15  and  16  limited  their  usefulness.  The  form  of  Equation  15  is  such 
that,  to  employ  the  un  du  integration,  the  exponent  m  must  be  less  than 
-2.  This  is  usually  the  case  in  bare-charge  craters,  where  the  well-defined 
lips  make  for  a  fairly  steep  decay  of  the  ejecta  field.  For  the  projectile 
craters,  however,  a  lip  is  almost  nonexistent,  and  fewer  than  half  of  the 
events  on  which  sampling  was  done  resulted  in  equations  that  could  be  inte¬ 
grated  by  means  of  Equation  16.  For  such  cases,  a  second  approach  to  ejecta 
areal  density  has  been  used.  This  equation  states 

6  =  A  e-BR  (17) 

where  A  and  B  are  constants.  Experience  with  large  bare-charge  craters 
has  shown  that  Equation  17  underestimates  the  ejecta  profile  in  the  crater 
lip,  but  it  was  not  known  whether  or  not  this  would  be  a  problem  with  the 
munitions  craters.  In  view  of  the  need  for  additional  data,  those  events  on 
which  ejecta  sampling  was  conducted  but  which  could  not  be  evaluated  by  means 
of  Equations  15  and  16  were  placed  in  the  form  of  Equation  17.  For  these 
events,  integration  by  parts  was  necessary  to  find  ejecta  weight.  Thus, 

We  «  2tt  I  A  e  R  dR  (18) 

a 

96.  Table  5  shows  the  are.il  density  equations  applicable  to  the  various 
TTC  events.  In  most  cases,  evaluation  by  Equations  17  and  18  resulted  in 
smaller  ejecta  quantities  than  for  Equations  15  and  16.  Based  upon  volumetric 
analyses,  it  appears  that  Equations  15  and  16  give  the  better  estimates;  how¬ 
ever,  all  data  have  been  included  in  the  analysis  leading  to  the  selection  of 
a  ve/va  fraction  for  the  TTC  shots — a  fraction  that  was  applied  to  all 
static-fire  events  for  which  ,no  better  information  was  available. 

97.  The  TTC  study  also  included  explosion-cloud  (dust  plus  explosion 
gases)  data.  While  there  is  no  attempt  here  to  address  cloud  dynamics. 

Table  5  does  include  maximum  recorded  areas  of  cloud  obscuration  measured  fr^m 
ground  stations  at  some  time  (seconds)  after  detonation.  Also  shown  are 
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vegetation  and  canopy-cover  descriptions,  since  these  may  influence  both  the 
cloud  and  ejecta  deposition. 

98.  One  further  comment  on  the  TTC  data  is  in  order:  placement  of  the 
projectiles  was  accomplished  by  means  of  prefabricated,  wooden  stands  to 
provide  a  firm  base  and  the  desired  angle.  It  seems  very  likely  that  the 
stands  absorbed  some  of  the  fragment  kinetic  energy  and  thus  probably  modified 
crater  formation  and  ejecta  deposition  to  an  unknown  degree  (but  believed  to 
be  small) . 

Uncased  charges 

99.  Bare-charge  craters  intended  to  simulate  projectile  detonations 
provide  the  final  source  of  ejecta  data.  As  with  projectiles,  the  preferred 
method  of  ejecta-volume  calculation  is  by  ejecta  sampling,  available  in  a  num¬ 
ber  of  the  bare-charge  tests.  Where  this  method  is  used,  areal-density 
Equation  15  is  shown  in  Tables  1-4  and  6-8,  along  with  calculated  ejecta 
weight  (Equation  16)  and  V  ,  found  by  dividing  W  by  wet  density.  It 
should  be  noted  here  that  questionable  results,  resulting  from  anomalous  sam¬ 
ples  have  been  discarded. 

100.  The  second  method  of  ejecta-volume  calculation  is  by  volumetric 
balance  (Equation  14),  discussed  in  paragraph  87;  this  method  is  possible 
where  true-crater  measurements  have  been  made.  In  the  case  of  bare-charge 
craters,  however,  a  better  understanding  of  crater  volumes  enabled  some 
approximations  that  could  not  be  reliably  made  for  projectile  craters — 
approximations  which  permitted  Vg  estimates  that  would  not  otherwise  have 
been  possible.  These  approximations  began  with  a  rearrangement  of 
Equation  14: 


Vfb  'Vc  +  Vfl> 


(19) 


The  terra  (V  -  V  )  can  be  substituted  for  Vr,  ,  since  fallback  comprises 
t  a  rb 

the  difference  between  the  true  and  apparent  crater  profiles.  It  can  be 
reduced  to  its  in  situ  volume  by  using  a  bulking  factor  of  1.5,  as  an  approx¬ 
imation.  A  small-scale  test  of  surface  detonations  (Chew,  Rooke,  and  Pitman 
1968)  actually  measured  both  V  and  V  ,  the  only  such  measurement  known 

C  X  Jt 

to  the  authors.  For  dry  sand,  it  was  found  that  V  <  0.02  V  ,  but  that  (V 

c  t  c 


+  vft>  -  0.27  Vt 


While  the  mechanisms  of  compaction  and  flowage  are 
expected  to  vary  considerably  with  differences  in  charge  geometry  and  cratered 


medium,  their  relation  to  true  crater  volume  should  be  fairly  constant  for  the 
ranges  of  conditions  in  this  study;  therefore,  this  assumption,  along  with  the 
foregoing  substitution,  permits  the  following  form  of  Equation  7: 


V 

e 


V 


0.06  V  +  0.67  V 
t  a 


0.27  V 


(20) 


Also, 


V  =  0.73  V„  -  0.67  V,,  (21) 

e  t  fb 

Both  Equations  20  and  21  have  been  used  in  the  data  tables. 

101.  The  foregoing  approaches  have  been  used  to  calculate  or  estimate 

V  and  V  /V  wherever  possible.  These  and  V  /V  relations  that  could  be 
e  6  d.  e  a 

found  in  the  bibliographic  literature  then  became  the  basis  for  estimation  of 

ejecta  quantities  where  only  apparent  crater  volume  was  known.  As  a  rule,  the 

search  was  restricted  to  the  charge  size  limitations  of  this  study,  since 

there  are  some  indications  that  these  relations  change  with  increasing  charge 

sizes.  Where  they  could  be  found,  measured  relations  were  recorded  for  each 

test  site.  As  would  be  expected,  however,  these  data  were  often  incomplete; 

and  in  some  cases,  there  was  lack  of  agreement  between  apparently  similar 

soils  at  different  test  sites,  so  that  the  final  selection  of  the  V  /V 

e  a 

relation  sometimes  became  a  matter  of  judgment.  Table  14  lists  the  relations 
derived  from  the  literature  search;  these  are  discussed  in  general  terms  in 
the  following  paragraphs.  Note  that  the  table  includes  some  small-scale 
studies  in  which  charge  weights  are  on  the  order  of  a  fraction  of  a  gram. 

These  are  viewed  with  somewhat  the  same  suspicion  as  large-scale  shots,  where 
charges  may  be  measured  in  tons.  These  small-scale  studies  were  very  care¬ 
fully  done,  however,  and  the  resulting  relations  are  useful  in  filling  the 
gaps  in  this  study. 

102.  Note,  in  Table  14,  that  the  silt/sand  and  clay/sand  combinations 

at  the  several  test  sites  differ  considerably  in  the  V  /V  relation.  While 

e  a 

relations  at  the  two  Fort  Carson  sites  were  quite  consistent  for  a  range  of 
charge  geometries,  the  tailored  soils  in  the  Fort  Polk  experiments  varied 
rather  widely.  The  WSMR  dry  silty  sand  showed  even  more  variation,  looking 


60 


more  like  the  dry  sand  at  the  Eglin  Air  Force  Base  site.  The  WSMR  data  repre¬ 
sent  an  average  of  three  shots  on  which  ejecta  sampling  was  done;  two  such 
shots  were  discarded  because  the  relations  were  obviously  invalid.  Although 
not  entirely  credible,  the  dry  sand  and  dry  silty  sand  data  seem  to  show  that 
a  high  fraction  does  indeed  exist  for  these  conditions.  It  is  specu¬ 

lated  that  this  may  be  caused  in  part  by  elastic  rebound,  in  which  the  bottom 
of  the  true  crater  "bounces"  upward  after  passage  of  the  compressional  shock 
wave,  leaving  a  true-crater  mound  in  which  in  situ  density  has  been  reduced  to 
the  point  that  overall  "compaction"  may  actually  be  negative.  This  event 
obviously  reduces  crater  size  and  may  lead  to  the  unusually  high  V  /Vfl 
fraction. 

103.  It  might  be  thought  that  an  increase  in  moisture  would  reduce  the 

V  /V  relation  in  soil,  since  it  would  be  expected  to  increase  both  the  plas- 
e  a 

ticity  and  cohesion  of  the  cratered  medium;  the  former  would  be  expected  to 

result  in  a  larger  crater  and  the  latter  in  a  reduction  in  ejected  material. 

The  20-ton  FLAT  TOP  II  and  III  shots  seemed  to  show  this.  FLAT  TOP  III  was 

identical  to  its  predecessor  except  for  soil  moisture,  which  had  approximately 

doubled,  resulting  in  a  crater  that  was  larger  by  more  than  50  percent,  but  in 

which  the  V  /V  relation  was  reduced  nearly  40  percent  (Ahlers  and  Miller 
e  a 

1966,  Rooke  and  Davis  1966).  While  this  holds  true  for  some  of  the  entries  in 
Table  14,  however,  it  is  not  seen  consistently. 

104.  Finally,  it  should  be  noted  that  in  some  cases  in  the  data  tables 
(the  BETS  experiments),  ejecta  relations  are  based  upon  shot-to-shot  moisture 
conditions  and  thus  may  seem  to  disagree  with  the  relations  in  Table  14.  In 
other  cases,  as  for  wet  sand,  the  relations  are  the  result  of  graphical  inter¬ 
polation  of  charge  position.  Table  14  gives  recommended  vg/va  fractions 
(after  rejection  of  invalid  measurements)  and  the  basic  statistics  associated 
with  each  entry,  where  appropriate. 

105.  Figure  24  graphs  V  /V  versus  R  for  the  bare-charge  data  of 
Table  14,  All  data  are  for  sands  and  fine-grain  soils,  the  main  difference 
between  cratered  media  being  that  of  moisture  content.  The  data  scatter, 
which  is  especially  prevalent  in  the  drier  soils,  seems  to  preclude  a  more 
elaborate  analysis.  Based  upon  a  visual  fit,  there  appears  to  be  a  slight 
downward  trend  as  the  charge  position  is  moved  upward.  For  the  surface  to 
surface-tangent  regime,  the  ratio  V  /V  lies  between  0.75  and  0.80.  Judg- 
ment  was  used  in  applying  these  ratios  and  those  of  Table  14  to  events  for 
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Figure  24.  Graphical  representation  of  ejecta  volume/ 
apparent  crater  volume  ratios  of  Table  14 


♦40 


which  no  ejecta  measurements  were  available.  In  some  cases,  averages  from 

different  sites  were  used.  Where  used,  selected  ratios  are  shown  under  the 

"Ejec  Cal"  column. 

Comparison  of  ejecta  from 
cased  and  uncased  charges 

106.  In  Table  13,  all  projectile  craters — static  and  live  fire — were 
examined,  and  the  resulting  Ve/Vg  =  0.52  was  used  to  estimate  ejecta  quan¬ 
tities  for  the  projectile  craters  of  Table  9  except  Shot  T-27,  for  which  a 
volumetric  balance  was  available.  However,  Table  5  was  not  changed  from  its 
listing  in  the  Phase  I  report  where  the  Vg/Va  ratio  equaled  0.53,  a  differ¬ 
ence  judged  not  to  be  significant  in  view  of  the  uncertainty  of  results 
obtained  in  the  Table  13  sample.  This  treatment  minimizes  the  Shot  T-27 

results,  which  are  considered  low.  Comparison  of  V  /V  ratios  for  cased  and 

e  a 

uncased  charges  indicates  that  the  ratio  for  the  former  is  only  about  two- 
thirds  that  of  the  latter. 
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Bare-Charge  Simulation  of  Munitions 


107.  Although  not  one  of  the  objectives  of  this  study,  it  might  be 

appropriate  to  discuss  some  of  the  findings  that  may  prove  helpful  to 

researchers  attempting  to  simulate  cratering  effects  of  munitions  using 

uncased  (bare)  charges  or  attempting  to  correlate  effects  observed  in  past 

experiments.  An  example  may  serve  to  illustrate  this  point. 

a.  Suppose  the  DOT  I  study  (Table  7)  is  reexamined  to  see  how 
closely  the  B-series  surface-tangent  (ST)  shots  simulated  the 
crater/ejecta  effects  of  a  155-mm  projectile  fired  in  an  ST 
geometry  approximating  live  fire.  Note  that,  in  terms  of  TNT 
equivalence,  charge  weights  are  quite  close:  6.46  kg  for  the 
bare  charges  versus  6.62  kg  for  the  155  mm.  Va  from  Rows 

404-413  in  Table  7  may,  of  course,  be  averaged,  or  an  estimate 
from  Figure  5  and  its  accompanying  equation  may  be  obtained: 

V  *  0.07  W0,78  (22) 

a 

Choosing  the  latter  course  tends  to  smooth  out  anomalous  mea¬ 
surements  such  as  the  suspiciously  large  V  of  Rows  407.  If 

cl 

it  is  decided  to  solve  the  equation,  the  following  is 
obtained: 


Va(bc)  '  °-07  «-«)°-78 

=  0.30  m3  (23) 

where  '7a(^c)  is  the  apparent  crater  volume  for  the  bare 

charge.  This  could  have  been  estimated  directly  from  the 
graph.  If  V  's  had  been  averaged,  0.34  m  would  have  been 

3 

obtained.  Following  a  similar  procedure  with  the  data  avail¬ 
able  on  the  155-mm  shell,  the  measurements  in  Tables  5  and  9 
may  be  averaged,  or  any  or  all  of  Figures  22  and  23  and 
Table  11  can  be  referred  to.  If  it  is  decided  to  use  the 
"live  plus  selected  static  firings"  and  the  Y  =  AX  fit, 
then 


V 


a(mun) 


=  0.071W1*07 
=  0.54  m3 


(24) 


where  ^a(mun)  refers  to  apparent  crater  volume  for  the  muni¬ 
tion.  Regarding  ejecta,  note  that  the  best  evidence  now 

available  indicates  that  the  V  /V  ratio  for  projectile  era- 

e  a 

ters  is  significantly  smaller  than  for  bare  charges.  Table  14 
or  Figure  24  can  be  consulted  for  a  bare-charge  ratio  and 
Table  13  for  projectiles.  Unfortunately,  Table  14  does  not 
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offer  any  clear-cut  choices  In  this  case;  therefore.  Figure  24 

Is  resorted  to,  which  gives  an  estimated  bare-charge 

V  /V  *  0.78  .  Table  13  gives  the  projectile  V  /V  =  0.52  . 
e  a  e  a 

Thus,  without  considering  dust  lofted  by  fragments  impacting 

beyond  the  crater, 


and 


V  .  -  0.78  V  . 
e(bc)  a(bc) 


0.23  m" 


V  ,  .  =  0.52  V  .  . 

e(mun)  a(mun) 


=  0.28  m' 


(25) 


(26) 


All  things  considered,  it  appears  that  a  155-mm  shell  fired  in 
a  fuze-quick  configuration  will  produce  a  slightly  larger 
quantity  of  ejecta  than  a  DOT  I,  B-series  ST  charge. 

b.  Continuing  with  the  example,  if  it  is  decided  to  design  a  fol¬ 
low-on  experiment  that  will  duplicate  the  DOT  I  B-series  ST 
shots  except  for  charge  weight  and  will  simulate  the  155-mm 
conditions  discussed  above  as  closely  as  possible,  it  can  be 
reasoned  that  more  ejecta  need  to  be  produced  from  the  bare- 
charge  crater  by  a  factor  of  0.28/0.23  *  1.22  .  Therefore, 
for  the  new  bare  charge  Nc.  2(bc2), 


To  find 
Va(bc2) 


1.22  Ve(bc)  -  0.78  Va(bc2) 

1.22  (0.23)  =  0.78  Va(bc2) 

V  ..  =  0.36  m3 
a(bc2) 

the  charge  weight  that  will  provide  the  required 
use  the  graph  of  Figure  6  or  solve 

t 

0.07  W0,78  =  0.36 


(27) 


W  =  8.13  kg  (TNT  equivalent)  (28) 

For  C-4,  with  a  95-percent  TNT  equivalence  in  cratering, 
actual  charge  weight  should  be  (1/0.95) (8. 13)  =  8.54  kg  . 

108.  In  the  above  example,  the  munitions  crater  is  larger  than  the 
bare-charge  crater,  a  result  that  is  consistent  with  Figure  23  in  the  charge- 
weight  range  above  5  kg.  However,  it  must  be  borne  in  mind  that  the  munitions 
data  base  is  small  compared  with  that  for  bare  charges  and  that  a  larger  num¬ 
ber  of  munitions  observations  under  more  varied  conditions  might  lead  to  a 
different  result.  Thus,  the  significance  of  the  differences  in  the  example  is 
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open  to  question.  The  reader  is  cautioned  that  the  comparative  procedure 
described  above  is  primarily  for  illustrative  purposes  and  should  not  be  con¬ 
strued  to  imply  that  the  differences  are  firmly  established.  Additional  data 
particularly  additional  artillery  data,  are  needed  for  this  purpose. 


PART  V:  SUMMARY  AND  CONCLUSIONS 


Crater  Shape  Factors 

109.  Shape  factors  have  been  established  for  both  bare-charge  and  muni¬ 
tions  craters  for  the  calculation  of  crater  volumes.  While  the  former  shows 
little  deviation  and  can  be  used  with  confidence  to  calculate  crater  volumes, 
the  latter  varies  widely  with  type  of  munition  and  test  conditions.  However, 
the  mean  value  of  the  munitions  shape  factors  is  very  close  to  that  for  bare 
charges,  so  that  bare-charge  factors  may  be  used  for  munitions  craters  when 
pattern  information  is  unavailable.  These  factors  are  0.45  and  0.40  for 
apparent  and  true  craters,  respectively.  Combining  shape  factors  with  it 
terms  in  the  equation  gives 


V  =  1.41  r2d  (29) 

a  a  a 

and 

Va  -  1.26  r2dt  (30) 


VBSSFs 


Apparent  craters 

110.  A  detailed  analysis  of  the  data  base  contained  in  Tables  1-9  shows 
that,  for  a  general  equation  relating  crater  volume  V  to  charge  weight  W  , 

V  «  Wn  (31) 

where  n  is  an  exponent.  The  value  of  n  should  be  0.87  rather  than  the 
currently  used  1.111.  The  actual  apparent  crater  VBSSFs  for  various  charge 
positions  and  soil  types  may  be  read  from  the  graphs  shown  in  Figures  10-13. 
True  craters 

111.  Regarding  true  craters,  the  exponent  on  charge  weight  is  unity  or 
1.0.  The  true  crater  VBSSFs  can  be  read  from  Figure  19  for  various  charge 
positions.  In  Figure  19,  there  is  a  tendency,  though  not  well  defined,  for 
clay  data,  particularly  moist-to-wet  clay,  to  plot  in  the  upper  portion  of  the 
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data  envelope,  followed  by  silts,  sands,  and  rocks,  the  final  plotting  near 
the  bottom  of  the  data  envelope. 


Ratio  of  True  and  Apparent  Crater  Volumes 


112.  The  ratio  of  the  true  crater  volume  to  the  apparent  crater  volume 

(V  /V  )  for  various  charge  positions,  without  regard  to  soil  type,  is  approxi- 
C  3 

mated  by: 


Charge  Position 

Above  or  Below  (-) 

Ground  Surface,  R 
_ c 

£1.0 

0 

-1.0 

-2.0 


V£/V 


a 


1.0 

1.0  -  1.3 
1.3  -  1.7 


1.0  -  3.0 


Effect  of  Vegetation  on  Crater  Volume 


113.  Craters  formed  by  explosions  on  bare  ground  are  usually  about 
10  percent  larger  bv  volume  than  craters  formed  in  a  similar  soil  but  with  a 
grass  cover  that  is  less  than  30  cm  high.  In  such  a  case,  it  is  presumed  that 
coverage  will  range  from  50  to  100  percent.  Where  the  vegetative  coverage  is 
nearly  100  percent  and  consists  of  taller  grasses,  bushes,  and  undergrowth,  as 
well  as  small  trees,  the  crater  in  the  vegetated  areas  will  be  70  to  80  per¬ 
cent  the  size  (by  volume)  of  a  crater  over  similar  soil  but  without 
vegetation. 


Effect  of  Soil  Moisture  on  Crater  Volume 


114.  Within  the  limitation  of  the  data,  it  appears  that  soils  that  are 
wet,  but  not  saturated  to  the  point  that  slope  failures  occur  around  the 
crater  wall,  generally  produce  craters  that  are  somewhat  larger  than  dry  or 
slightly  moist  soils.  For  wet  sand,  the  linear  dimensions  of  a  crater  are 
about  5  percent  larger  than  for  dry  or  slightly  moist  sand.  For  silts  and 
clays,  the  linear  dimensions  of  craters  are  typically  about  10  percent  larger 


in  wet  soils  than  in  those  that  are  dry  or  slightly  moist.  These  increases 
translate  into  a  15-percent  volume  increase  for  sand  and  a  25-  to  30-percent 
increase  for  silts  and  clays  when  comparing  dry  or  slightly  moist  with  wet 
conditions.  As  to  rock,  no  effect  could  be  discerned  for  varying  moisture 
content  within  the  yield  range  of  this  study. 

Cratering  the  HE  Projectiles 

115.  Craters  formed  by  HE  projectiles  detonated  in  live  fire  or  in  a 
geometry  to  simulate  live  fire  differ  markedly  from  comparable  bare-charge 
craters  in  shape,  but  not  in  volume.  As  to  comparisons  between  live  and 
static  fire  of  such  projectiles,  there  appears  to  be  little  difference  in 
overall  crater  size  and  shape. 


Ejecta  Assessment 

116.  Tables  1-9  contain  55  events  in  which  there  are  sufficient  data  to 
permit  direct  calculation  (by  integration  of  areal-density  measurements)  of 
|  crater-ejecta  volumes.  This  number,  representing  about  10.6  percent  of  the 

data  base,  is  approximately  evenly  divided  between  bare-charge  and  munitions 
events.  Another  16.5  percent  of  the  data  base  is  suitable  for  inference  of 
ejecta  quantities  by  means  of  volumetric  balancing  of  the  crater,  all  such 
I  balancing  being  based  upon  certain  assumptions.  Thus,  in  all,  about  27  per¬ 

cent  of  the  total  data  base  has  been  used  to  derive  ejecta/apparent-crater 
volume  ratios  (V^V^)  with  which  to  estimate  ejecta  volumes  for  the  remaining 
events  on  a  shot-by-shot  basis. 

|  117.  For  bare  charges,  24  V  /V  ratios  were  derived,  representing 

I  e  a 

most  of  the  soil  conditions  encountered  in  the  data  base.  Charge  positions 
were  predominantly  surface  to  surface-tangent  (0  ^  Rc  £  1.0).  While  the  data 
scatter  seems  to  preclude  rigorous  analysis,  its  general  trend  indicates  that 
j  0.75  <  V  /V  <  0.80  in  the  near-surface  regime,  with  a  slight  downward  trend 

as  charge  position  is  moved  from  below-surf ace  to  above-surface.  Thus,  in 
bare-charge  simulation  of  munitions,  roughly  three-fourths  of  the  apparent 
crater  volume  is,  on  the  average,  lofted  as  ejecta. 

118.  In  the  case  of  munitions,  only  the  105-  and  155-mm  artillery  pro¬ 
jectile  craters  were  usable  for  calculations  of  V  /V  ,  and  these  for  only 

6  3 

f 

I 


three  distinctly  different  soil  conditions.  Of  these,  a  single  live-fire 
event  was  available  for  the  vg/va  calculation.  Thus,  in  estimating  crater- 
ejecta  quantities  for  munitions,  a  limited  number  of  charge  and  soil  con¬ 
ditions  necessarily  served  as  the  basis  for  the  remaining  events.  The 

munitions  ratios  were  lower  than  for  bare  charges,  with  V  /V  -  0.52  ,  or 

e  a 

about  one  standard  deviation  below  the  inspected  mean  shown  in  Figure  24. 

However,  to  this  ratio  should  be  added  whatever  surface  material  is  dis- 

« 

sociated  and  lofted  by  shell  fragments  impacting  outside  the  crater;  this 
quantity  is  not  now  known. 

Bare-Charge  Simulation  of  Munitions 

119.  Using  the  established  data  base  and  the  crater/ejecta  relations 
derived  therefrom,  a  rational  approach  toward  simulation  of  munitions  with 
bare  charges  is  demonstrated  in  paragraph  107.  The  illustration  tends  to  sub¬ 
stantiate  field  observations  that  bare  charges  produce  smaller  craters  than  do 
comparable  munitions,  but  also  shows  that,  for  ejecta/dust  production,  this 

difference  is  narrowed  by  the  higher  V  /V  ratio  exhibited  by  bare  charges. 

0  0 
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PART  VI :  RECOMMENDATIONS 


120.  In  performing  the  analysis  and  in  formulating  the  conclusions  that 
were  derived  from  the  data  base,  one  of  the  more  glaring  difficulties  was  the 
continuing  problem  of  data  scatter.  Sources  of  data  scatter,  in  the  main,  are 
believed  to  be:  inconsistencies  in  crater  measuring  techniques,  i.e.,  the 
definition  of  the  crater  profile  (particularly  the  true  crater  profile);  the 
inherent  or  subtle  differences  in  reportedly  similar  soils  at  different  sites, 
including  differences  in  soil  composition,  moisture  content,  degree  of  satura¬ 
tion,  etc.;  differences  in  conversion  factors  regarding  explosion  yield  among 
different  types  of  explosives;  and  finally,  the  need  to  group  similar  shot- 
geometry  ranges,  since  data  for  unique  shot  geometries  were  not  available  in 
sufficient  quantity  to  permit  confident  analysis. 

121.  Wherever  the  degree  of  data  scatter  is  found  to  be  unacceptable, 
then  additional  testing — carefully  designed  and  controlled — is  recommended. 
However,  the  experimenter  is  advised  that  the  most  carefully  controlled  cra¬ 
tering  experiments  have  shown  a  10-  to  15-percent  scatter  in  linear  dimen¬ 
sions,  which  translates  to  roughly  30-percent  scatter  in  volume  for  a  specific 
soil  type,  a  specific  charge  position,  and  a  single  charge  weight. 

122.  Additional,  carefully  controlled  experiments  are  needed  to  better 
define  quantities  and  distribution  parameters  relative  to  crater  ejecta, 
especially  those  pertaining  to  both  live  and  statically  fired  munitions  and 
their  comparisons  with  bare  charges.  Where  possible,  direct  measurements  of 
ejecta  are  preferred  for  this  purpose,  as  by  sampling,  since  this  provides 
both  ejecta  quantity  and  distribution.  The  sampling  techniques  employed  in 
the  experiments  reported  herein  are  generally  adequate,  but  would  probably 
provide  better  results  if  sampling  arrays  were  expanded.  This  technique  does 
not,  of  course,  apply  to  live  fire,  where  the  projectile  detonation  point  can¬ 
not  be  known  ahead  of  time;  here  only  volumetric  analysis  is  applicable. 
Nevertheless,  estimates  of  ejecta  volume  Vg  may  be  inferred  from  careful 
measurements  of  crater  profiles.  Using  the  volumetric  relations  discussed  in 
paragraph  87  and  assuming  the  volume  of  compaction  that  surrounds  the  true 
crater  to  be  small,  then 


V  = 
e 


(32) 
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where  V  is  the  ejecta  volume,  V  is  the  apparent  crater  volume,  and 

e  a  f  £ 

is  the  volume  of  plastic  flowage,  which  is  the  same  as  the  volume  of 
upthrusted  (in  situ)  material  in  the  crater  lip,  as  shown  in  Figures  2  and  3. 

123.  In  considering  dust  loadings  resulting  from  artillery  fire,  both 
technical  and  tactical  realities  must  be  kept  in  mind.  Shell  fragments  cause 
dust  lofting  outside  fuze-quick  craters  in  dry  soils,  as  well  as  ejecta  from 
inside  the  craters.  Therefore,  crater-ejecta  volume  alone  will  not  account 
for  the  total  amount  of  dust  in  such  cases.  In  addition,  airbursts  are  fre¬ 
quently  used  against  soft  targets,  such  as  personnel,  and  these  bursts  may 
result  in  very  significant  dust  loadings  without  the  formation  of  any  crater. 
Here  again,  additional  field  work  is  probably  required  to  obtain  realistic 
input  to  models.  It  may,  however,  be  possible  to  obtain  estimates  of  dust 
loadings  from  artillery  fragments  impacting  the  ground  surface  by  reviewing 
and  analyzing  existing  motion-picture  photography  of  artillery  fire  that  is 
available  at  the  US  Army  Field  Artillery  School,  Fort  Sill,  Okla.  Prior  to 
incurring  the  expense  of  additional  field  work,  it  is  recommended  that  this 
approach  at  least  be  investigated. 

124.  Lastly,  cratering  effects  of  conventional  aerial  bombs  also  need 
additional  study.  This  is  another  munition  that  is  commonly  used  on  the  bat¬ 
tlefield.  General-purpose  bombs  are  relatively  thinly  cased  and  employ  dif¬ 
ferent  fuzing  from  artillery  projectiles;  hence,  penetration  and  blast  effects 
differ,  also.  Further,  armor-piercing  and  semiarmor-piercing  bombs  have  their 
own,  separate  characteristics  that  influence  cratering.  Some  bomb-crater  data 
are  available  from  past  studies,  but  their  suitability  for  ejecta  assessment 
is  questionable.  In  this  study,  the  upper  limit  of  energy  yield  is  sufficient 
to  include  bombs,  but,  for  the  reasons  cited  above,  crater/ejecta  predictions 
derived  herein  may  not  represent  well  the  actual  results  for  these  munitions. 
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KEY  TO  VEGETATION/CANOPY  DESCRIPTIONS,  TABLES  1-9 

AZ  -  Short  grasses,  blue  grama,  Russian  thistle,  lamb's  quarters,  and/or 
prickly  pear 

BL  -  Black  sand 

BS  -  Bare  soil 

C  -  Covered 

CG  -  Cut  grass 

DS  -  Dry  sand  (white) 

G  -  Gynerium  Sagittatvm 

GH&F  -  Grasses,  herbs,  and  forbs 

HG  -  High  grass 

LC  -  Low  canopy 

LG  -  Low  grass 

M  -  Morning  glory 

0  -  Other 

P  -  Panicum  sp  (1  to  2  m  tall) 

S  -  Spider  lily 
U  -  Uncovered 
WS  -  Wet  sand 

(50)  -  Indicates  percent  of  ground  covered  (in  this  case  50  percent) 


KEY  TO  COLUMN  HEADINGS  AND  ENTRIES,  TABLES  1-9 


1 

Row 

Defines  row  number 

2 

REF/SR 

Ref erence/subref erence 

3 

TEST  P/S 

Test  program/series 

4 

A  OBS 

Average  of  observations 

5 

SHOT  ID 

Shot  identification (s) 

6 

X  COORD 

X-coordinate 

7 

Y  COORD 

Y-coordinate 

8 

SOIL  DES 

Soil  description 

9 

uses 

Unified  Soil  Classification  System 
designation 

10 

SA,  1 

Percentage  of  sand  in  soil 

11 

SI,  % 

Percentage  of  silt  in  soil 

12 

CL,  % 

Percentage  of  clay  in  soil 

13 

MC,  Z 

Moisture  content,  percent 

14 

DESA,  % 

Degree  of  saturation,  percent 

15 

TYPE  CH 

Type  explosive  (charge) 

16 

CH  WT,  kg 

Charge  weight,  W  ,  kilograms 

17 

F  TNT,  kg 

Equivalent  charge  weight,  TNT,  kilograms 

18 

TM  C/M 

Type  munition,  caliber /model 

19 

MO  D/A 

Munition  orientation,  description/angle 
Notation  used  in  column  entries: 

S  -  Surface  (half-buried)  geometry 
ST  -  Surface-tangent 
STB  -  Surface-tangent-below,  or 

buried  one  projectile  length 
1/3BUR90  -  Projectile  buried  one-third  of 
its  length  and  at  90-deg 
angle  from  horizontal 
1/3BUR45  -  Projectile  buried  one-third  of 
its  length  and  at  45-deg  angle 
from  horizontal 

20  H/DOB,  R  Charge /munit i on  height  or  depth  of 

burst  in  equivalent  TNT  charge  radii 


KEY  TO  COLUMN  HEADINGS  AND  ENTRIES,  TABLES  1-9 
Continuation  Sheet  of  Each  Table 


Column  No. 

Heading 

Meaning 

1 

Row 

Defines  row  number 

2 

Shot  ID 

Shot  identif ication(s) 

3 

r  ,  m 
a 

Apparent  crater  radius,  metres 

4 

d  ,  m 
a 

Apparent  crater  depth,  metres 

5 

V  ,  m3 
a 

Apparent  crater  volume,  cubic  metres 

6 

r  ,  m 

True  crater  radius,  metres 

7 

dt  *  ”, 

True  crater  depth,  metres 

8 

3 

Vt  »  m 

True  crater  volume,  cubic  metres 

9 

*  3 

V  m  /W 
a  s 

Apparent  crater  volume /charge 

n  87 

weight  to  the  0.87  power  (V  /W  ) 

a 

10 

V  /W 
a  n 

Apparent  crater  volume /charge 

weight  to  the  1.111  power  (V  /W*‘^') 

11 

Vt  m3/kg 

True  crater  volume /charge  weight, 

m3/kg 

12 

V  /Wn 

True  crater  volume /charge  weight  to  the 
1.111  power  (m3/kg*'*^) 

13 

V  ,  m3 
e 

Volume  of  ejecta,  cubic  metres 

14 

We  .  kg 

Weight  of  ejecta,  kilograms 

15 

Ejec  Cal 

Basis  of  ejecta  volume  calculations — 
equation  based  on  sampling,  volumetric 
balance,  or  e jecta-apparent  crater 
volume  ratio 

16 

D  C  OBS 

2 

Dust  cloud  observation,  m  /sec.  Maximum 
target  area  obscured  at  t  seconds 
after  detonation. 

17 

K 

s 

Crater  shape  factor  (apparent  or  true 
crater) 

18 

V/C  Des 

Vegetation/canopy  description 
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KEY  TO  SOIL  DESCRIPTIONS,  TABLES  1-9 


Soil  Description 

Nomenclature 

Meaning 

DDALLUV 

Dry  desert  alluvium 

DTMSand 

Dry-to-moist  sand 

IC1 

Inorganic  clay 

ICl-SaM 

Inorganic  clay-sand  mixture 

ISiCISi 

Inorganic  silt — clayey  silt 

Kaolin 

Fat  clay 

M-Loess 

Moist  loess 

MSaSilt 

Moist  sandy  silt 

MSaSiCl 

Moist  sandy  silty  clay 

MSiClay  (MSiCl) 

Moist  silty  clay 

MS i Sand 

Moist  silty  sand 

PGSa 

Poorly  graded  sand 

PUMSAND 

Pumice  sand 

Rock  B-G 

Rock,  basalt,  or  granite 

SANDST 

Sandstone 

SaCl 

Sandy  clay 

SaCl.P 

Sandy  clay,  plastic 

SaCISi 

Sandy  clayey  silt 

SaSi 

Sandy  silt 

SaSiCl 

Sandy  silty  clay 

SH-T-FG 

Shale,  tuff,  or  frozen  grou 

SiCl 

Silty  clay 

SiSand  (SiSa) 

Silty  sand 

SiSa  Mix 

Mixture  of  silt  and  sand 

UNISAND 

Uniform  sand 
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Cratering  Data  from  WES  Experiments  in  Wet  Sand  and  from 
Munitions/Bare-Charge  Equivalence  II  Study 
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Table  12 

Cratering  Results  for  IPS-  and  155-mm  HE  Projectiles  In  Live  or  Simulated  Live  Fire 
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Interval  calculations  based  upon  °n  -  1  for  95-percent  area  under  normal  curve  or  for  Student's  t-distributlon  for  5-percent  tall. 
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